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This review presents the elastic theory of low-dimensional (one- and two-dimensional) continua 
and its applications in bio- and nano-structures. 

First, the curve and surface theory, as the geometric representation of the low- dimensional con- 
tinua, is briefly described through Cartan moving frame method. The elastic theory of Kirchhoff 
rod, Helfrich rod, bending-soften rod, fluid membrane, and solid shell is revisited. The free energy 
density of the continua, is constructed on the basis of the symmetry argument. The fundamental 
equations can be derived from two kinds of viewpoints: the bottom-up and the top-down stand- 
points. In the former case, the force and moment balance equations are obtained from Newton's 
laws and then some constitute relations are complemented in terms of the free energy density. In 
the latter case, the fundamental equations are derived directly from the variation of the free en- 
ergy. Although the fundamental equations have different forms obtained from these two viewpoints, 
several examples reveal that they are, in fact, equivalent to each other. 

Secondly, the application and availability of the elastic theory of low-dimensional continua in 
bio-structures, including short DNA rings, lipid membranes, and cell membranes, are discussed. 
The kink stability of short DNA rings is addressed by using the theory of Kirchhoff rod, Helfrich 
rod, and bending-soften rod. The lipid membranes obey the theory of fluid membrane. The shape 
equation and the stability of closed lipid vesicles, the shape equation and boundary conditions of 
open lipid vesicles with free edges as well as vesicles with lipid domains, and the adhesions between 
a vesicle and a substrate or another vesicle are fully investigated. A cell membrane is simplified 
as a composite shell of lipid bilayer and membrane skeleton, which is a little similar to the solid 
shell. The equations to describe the in-plane strains and shapes of cell membranes are obtained. It 
is found that the membrane skeleton enhances highly the mechanical stability of cell membranes. 

Thirdly, the application and availability of the elastic theory of low-dimensional continua in nano- 
structures, including graphene and carbon nanotubes, are discussed. A revised Lenosky lattice model 
is proposed based on the local density approximation. Its continuum form up to the second order 
terms of curvatures and strains is the same as the free energy of 2D solid shells. The intrinsic 
roughening of graphene and several typical mechanical properties of carbon nanotubes are revisited 
and investigated based on this continuum form. It is possible to avoid introducing the controversial 
concepts, the Young's modulus and thickness of graphene and single-walled carbon nanotubes, with 
this continuum form. 
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I. INTRODUCTION 

We human beings live in a three-dimensional (3D) 
space which contains many geometric entities composed 
of atoms or molecules. The length scale of objects ob- 
served with our naked eyes is much larger than the dis- 
tance between nearest neighbor atoms or molecules in 
the objects. As a result, the objects can be regarded as 
continua. If one dimension of an object is much larger 
than the other two dimensions, such as a rod, we call it 
a one-dimensional (ID) entity. If one dimension of an 
object is much smaller than the other two dimensions, 
such as a thin film, we call it a two-dimensional (2D) en- 
tity. In this review, the term "low-dimensional continua" 
represents ID and 2D entities. 

Elasticity is a property of materials. It means that 
materials deform under external forces, but return 
to their original shapes when the forces are removed. 
Elastic theory, the study on the elasticity of contin- 
uum materials, has a long history-^ which records 



many geniuses such as Hooke (1635-1703), Bernoulli 
(1700-1782), Euler (1707-1783), Lagrange (1736-1813), 
Young (1773-1829), Poisson (1781-1840), Navier 
(1785-1836), Cauchy (1789-1857), Green (1793-1841), 
Lame (1795-1870), Saint- Venant (1797-1886), Stokes 
(1819-1903), Kirchhoff (1824-1887), and so on. Now 
elastic theory has been a mature branch of physics 
and summarized in several excellent textbooks! 2 -! 3 -^ 
Although the classical elastic theory is applied to 
macroscopic continuum materials, more and more facts 
reveal that it can be also available for bio- or nano- 
structures such as short DNA rings s&LS&l&lhl&lS. 
a-helical coiled coils ^ chiral filaments ,I5 1 i6 i r7 1 i8 1 i9 1 20 i 2i 
climbing plants ; 22 ! 23 bacterial flagella^ viral 
shells^* 2 ! bio-membranes^^i^ 3 .^^ 
zinc oxide nanoribbonS ) 37 i 38 ' 39 and carbon 
nanotubeS ) 40 i 41 ' 42 i 43 ' 44 i 45 ' 46 to some extent. 

This review presents the elastic theory of low- 
dimensional continua and its applications in bio- and 
nano-structures, which is organized as follows: In Sec.HIl 
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we briefly introduce the geometric representation and the 
elastic theory of low-dimensional continua including ID 
rod and 2D fluid membrane or solid shell. The free en- 
ergy density of the continua is constructed on the basis 
of the symmetry argument. The fundamental equations 
can be derived from the bottom-up and the top-down 
viewpoints. Although they have different forms obtained 
from these two standpoints, several examples reveal that 
they are, in fact, equivalent to each other. In Sec. lIIIl the 
application and availability of the elastic theory of low- 
dimensional continua in bio-structures, including short 
DNA rings, lipid membranes, and cell membranes, are 
discussed. We investigate the kink stability of short DNA 
rings, the elasticity of lipid membranes, and the adhe- 
sions between a vesicle and a substrate or another vesicle. 
A cell membrane is simplified as a composite shell of lipid 
bilayer and membrane skeleton. The membrane skeleton 
is shown to enhance highly the mechanical stability of 
cell membranes. In Sec. IIV1 the application and avail- 
ability of the elastic theory of low-dimensional continua 
in nano-structures, including graphene and carbon nan- 
otubes, are discussed. We propose a revised Lenosky lat- 
tice model and fit four parameters in this model through 
the local density approximation. We derive its contin- 
uum form up to the second order terms of curvatures and 
strains, which is the same as the free energy of 2D solid 
shells. The intrinsic roughening of graphene and several 
typical mechanical properties of carbon nanotubes are 
revisited and investigated by using this continuum form. 
Sec. [V] is a brief summary and prospect. 



the Frenet frame which satisfies the Frenet formula^! 

r r' = T, 
T' = K N, 

S N' = (-kT + tB), { ' 

{ B' = -tN, 

where the prime represents the derivative with respect to 
s. k and r are the curvature and torsion of the curve, 
respectively. 



C 




FIG. 1: Frenet frame {r;T,N,B}. 



The fundamental theory of curved tells us that the 
bending and twist properties of a smooth curve are 
uniquely determined by the Frenet formula JTJ). 

2. Surface theory 



II. FUNDAMENTALS OF GEOMETRIC AND 
ELASTIC THEORY ON LOW-DIMENSIONAL 
CONTINUA 

In this section, we describe the mathematical basis and 
the elastic theory of ID and 2D continua. 



A. Geometric representation of low-dimensional 
continua 

The ID continuum (rod) and 2D continuum (mem- 
brane or shell) can be expressed as a smooth curve and 
a smooth surface, respectively. 



1. Curve theory 

Fig. m depicts a curve C embedded in the 3D Euclid 
space. Each point in the curve can be expressed as a 
vector r and let s be the arc length parameter. At point 
r(s), one can take T, N, and B as the tangent, normal 
and binormal vectors, respectively, {r; T, N, B} is called 



Fig. [5] depicts a surface M embedded in the 3D Euclid 
space. Imagine that a mass point moves on the surface in 
the speed of unit and that a right-handed frame, which 
consists of three unit orthonormal vectors with two vec- 
tors always in the tangent plane of the surface, adheres 
to the mass point. Assume that the mass point is at posi- 
tion expressed as vector r and the frame superposes three 
unit orthonormal vectors {ei, 62,63} with 63 being the 
normal vector of surface M at some time s. When the 
mass point moves to another position r' at time s + As, 
the frame will superpose three unit orthonormal vectors 
{e^, e 2 , e 3 }. Thus we call the frame a moving frame and 
denote it as {r; ei, e 2 , e 3 }. 

If As — > 0, we define 

dr = lim (r' — r) = aiiei + 102^2* (2) 

and 

efe; = lim (e- - ej) = Wwe,-, (i = 1, 2, 3) (3) 

As— >0 

where u>\, u>2, and wy, (i, j — 1,2,3) are 1-forms, and '<f 
is the exterior differential operator . 48 ' 49 Here LO12 can be 
understood as the infinite rotation angle of vectors ei and 
e 2 around e3. Similarly, we can understand the physical 
meaning of the other w.y . It is easy to obtain Wy = — uji 
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FIG. 2: Moving frame {r;ei,e2,e3} of a surface M. 



from e; • ej = 8tj . Additionally, the structure equations 
of the surface can be expressed as i 48 ' 49 

du>i = CJ12 A 1V2, 

dui2 = u>2i A w\, (4) 
duij = u) ik A w kj (i, j = 1, 2, 3), 

and 




where 'A' represents the wedge production between two 
differential forms. The matrix ( ^ ^ ~\ is the represen- 
tation matrix of the curvature tensor 9\. Its trace and 
determinant are two invariants under the coordinate ro- 
tation around e3 which are denoted by 

2H = a + c and K = ac - b 2 . (6) 

They can be expressed as 2H = — (1/Ri + I/R2) and 
K = I/R1R2 by the two principal curvature radii R\ 
and i? 2 at each point. 

Consider a tangent vector m stemming from r. Let 
4> be the angle between m and ei. Then the geodesic 
curvature, the geodesic torsion, and the normal curvature 
along the direction of m can be expressed^ 

k g = (d(f) + V12) /ds, 

T g = bcos2(f> + (c — a) cos (/) sin 0, (7) 
k n = a cos 2 4> + 2b cos (f> sin <j> + c sin 2 (j), 

where ds is the arc length element along m. If m aligns 
with ei, then 4> — 0, k g — u^/ds, r g — b, and k n = a. 

B. Elastic theory of ID continua 

We will elucidate the elastic theory of rod with inex- 
tensible centerline. As shown in Fig. [3l let us simplify a 
rod as a curve r(s) with s being the arc-length parameter, 
and cut an infinitesimal element (shown in the magnified 
box) from the rod. There are forces and moments at the 



M+dM 




FIG. 3: Force and moment in ID rod. 



two ends of the element which originating from the in- 
teraction of other parts of the rod. F and M represent 
the force and moment vectors at point r(s), while F + c£F 
and M + c?M are the force and moment vectors at point 
r(s + ds). From Newton's laws, we can derive the force 
and moment balance equations: 

^F = 0^F' = 0, (8) 

and 

^M = 0^M' + TxF = 0, (9) 

where the prime represents the derivative with respect to 
s. One should add the constitutive relation and boundary 
conditions to make the above two equations closed. 

1. Kirchhoff rod theory 

A rod with rectangle cross section and centerline C is 
shown in Fig. [4] Take local coordinates {x\, X2, £3} with 
x\ and X2 paralleling respectively to the two edges of the 
rectangle, and X3 along the tangent of the centerline. N is 
the normal of curve C. Let {xi,X2,X3} denote the basis 
of the local coordinates and define n± = — x 2 • (dx 3 /ds), 
«2 = xi • (dx 3 /ds), and k 3 = x 2 • (dxi/ds). Viewed 
from geometrical point, k± and K2 describe the bending 
of the rod around axes x\ and X3, respectively, and K3 
represents the twist of the rod around axis 23. The free 
energy density G due to the bending and twist can be 
expressed as a function of Ki, K2, and K3. Expanding G 
up to the second order terms of K\, k 2 , and K3, we have 

G = 7+ ^( Kl -^) 2 + y (K2-S 2 ) 2 + ^( K 3-S 3 ) 2 , (10) 

where the constant 7 can be interpreted as the line ten- 
sion. R\ and R-2 are interpreted as the spontaneous 
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curvatures while R 3 the spontaneous torsion. Denote 
k = K1X1 + K2X2 + K3X3 and let be the angle between 
Xi and N. Then we have 



N = cos xi — sin x 2 , 
B = sin 6 xi + cos X2 , 



(11) 



where B is the binormal of curve C. From Eqs. fl} and 
(fTTj) . we can deriv o 50 ' 51 



k = «sin0xi + k cos 0X2 + (r + <t>') X3. 
Thus G can be also regarded as the function of n, r, i 



(12) 




FIG. 4: Rod with rectangle cross section. 



The moment vector is defined as^ 
dG dG 



_ 5G _ 

9k 9ki 1 



<9G 



a X l + a X 2 + a X 3; (13) 



dn 2 



dn 3 



which is called the constitutive relation. Eqs. (JSJ), © 
and (|13p with some boundary conditions form a group of 
closed equations. They are also available for the rod with 
cross section different from rectangle if only we take x\ 
and x 2 as the two principal axes of inertia. It should be 
noted that the equivalent form of these equations can be 
also obtained from the variational method. This method 
is called the top-down method while the former one via 
Newton's laws called the bottom-up method. 

The free energy of a rod with length L can be written 

as 



G(k, t; 0, 0') ds + Ti 



(14) 



where Thd comes from the contributions of two ends of 
the rod. The general Euler-Lagrange equations corre- 
sponding to Eq. (fT4")) are derived as 

G*-(GV)' = 0, (15) 
G" K + 2t(G;/k)' + G' t t'/k + (k 2 - t 2 )G k 

+2ktG t + Ktfi'Gp -kG = 0, (16) 
t'G k + 2tG' k - (kG t )' + {t 2 /k)G' t - (G;/k)" = 0, (17) 

where G^, Gp, G K and G T are the partial derivatives of 
G with respect to 0, cf>' , k and r, respectively. Addition- 
ally G' K = (G K )', G' r = (G r )', G" K = (G K )". The berief 



derivation of Eqs. lfT5|) - IT7)) is attached in Appendix O 
These equations have been employed to investigate heli- 
cal and twisted filaments. ■AS There might be a misprint in 
Eq. (7) of Ref. UH, corresponding to our above equation 
(fT7)) . because the dimension of its last term is different 
from that of other terms. 

Now we would give a typical example to reveal 
the equivalence relation between Eqs. ([H]),©, (| 13[) and 
Eqs. (fl"5|) - (fl~7f rather than prove it directly. Let us con- 
sider a rod with k\ = k 2 — fco, ^3 = 0) an d R\ = R 2 = 
R 3 = 0. The free energy density (flQ|) is simplified as 



G = (fco/2)(/s? + «?) + 7 - (k /2) K 2 + 7. 



(18) 



On the one hand, we have Mi = fco K i = fcoKsin0, M 2 = 
kaK 2 — kon cos 0, M 3 = from Eq. (j 13|) . The moment 
balance equation ((9J implies F\ = —koKiK 3 — koK 2 and 
F 2 = k^n'i — koK 2 K 3 . Substituting them into the force 
balance equation {HJ, we have F 3 = F30 — koK 2 /2 and 

k" - kt 2 + k 3 /2 - nF 30 /k = 0, (19) 
2tk' + kt' = 0, (20) 

where F30 is an integral constant which represents the 
line tension of the straight (k = 0) rod. On the other 
hand, we have G K . — k K, G$ — G& = G T = 0. Eq. (fT5|) 
is trivial while Eqs. (fT6|) and (fT7|) are, respectively, trans- 
formed into 



0, 



(21) 
(22) 



k" - kt 2 + k 3 /2 - jn/k 
2tk + kt' = 0. 



The above equations are the same as Eqs. (|19|) and (|20|) 
obtained from the force and moment balance conditions if 
only we take F30 = 7. Thus the equations obtained from 
the top-down and bottom-up methods are equivalent to 
each other. 

Substituting the free energy density (fT0|) into 
Eqs. (fT5| - (fT7|) . we obtain the so called shape equations 
of Kirchhoff rod as 



(fci — k 2 )n 2 sin 20 — 2k 3 {r 



2I 21 k = 0, (23) 



Ii(2k" + k 3 - 2kt 2 ) - 2 7 k + 2/12 {<j> 12 + t ) 
+2/ 2 i0" +2(ki- k 2 )[{4>'K sin 20)' + 0V sin 20] 
-Ik + 4A: 3 t[(t' + <(/')/ k]' + 2fc 3 (r' + 0")t'/« 
+k 3 K(T + (f>' - R 3 )(3t + 4>' + R 3 ) = 0, (24) 



I^t'k + 2tk') - k 3 [n(T + 0' - R 3 )]' 

-k 3 T 2 (T' + ( f > ")/ K -k 3 [(T' ■ + </>")/ K]" 

-2r[(fei - /c 2 )Ksin20 + J 2 i]0' - W = 0, 



(25) 



where I\ = ki sin + ^2 cos 0, I = kiRf + k 2 R 2 , 1\ 2 = 
k\K\ sin0 + k 2 R 2 cos 0, and /21 = ^2^2 sin0 — kiRi cos 0. 

We also suggest that gentle readers consult the work 
by Zhou et al^ s where the above equations (I23p -((25 |) 
and different kinds of boundary conditions are expressed 
in another representation with the aid of Euler angles. 
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2. Helfrich rod theory 

Hclfrich rod theory can be regarded as the fourth order 
Kirchhoff rod theory with circular cross section to some 
extent. The free energy density is expressed as^ 

G= ifc 2 K 2 + fc 3K 2 T+ifc 22 K 4 + ifc 4 (K' 2 + K 2 T 2 )+7, (26) 

where fc 2 , fc 3 , fc 22 and fc 4 are elastic constants while 7 is 
the line tension. It is noted that this free energy den- 
sity is the simplest stable form including the chirality 
term but without spontaneous curvature and torsion. It 
has been employed to investigate the circular DNA in 
Ref. H and the Euler-Lagrange equations corresponding 
to J Gds are given as: 

fc 2 (/« 3 /2 — kt 2 + k") — 7« 
+fc 3 (3K 3 r - 2kt 3 + 6k V + 2kt" + 6k" t) 
+fc 4 (5K 3 r 2 /2 - kt 4 + kk' 2 /2 - k 2 k" - k"" 
+6k"t 2 + 12k'tt' + 4ktt" + 3kt' 2 ) 
+fc 22 (3K 5 /4 - k 3 t 2 + 6kk' 2 + 3k 2 k") = 0, (27) 



closure of single- walled carbon nanotubes with 1,3- 
dicyclohexylcarbodiimide^ Fitting the experiment data 
with the worm-like chain model,— the persistence length 
is 800 nm for single-walled carbon nanotubes in the di- 
ameter of 1 nmj£ which is much smaller than the the- 
oretical value 33 /jm estimated in terms of the Young's 
modulus and thickness of single- walled carbon nanotubes 
in Ref. 0. 

(a) (b) 
M M 




FIG. 5: Bending-induced softening relation: (a) the first kind 
in expression of Eq. (|29[) ; (b) the second kind in expression of 

Eq. daij. 



fc 2 (2«:'r + kt') + fc 22 (K 3 r' + 6k 2 k't) 
+fc 3 (6K'r 2 + 6ktt' - 3/cV - 2k'") 
+fc 4 (4K'r 3 + 6ktV - 3/tVr - «V 
-Ak't" - 6k" V - 4k"'t - kt'") = 0. (28) 

Here we will not go on the more higher order Helfrich 
rod theory, on which gentle readers can consult Refs. I20I 
andU 

3. Theory of bending-sofien Rod 

There are two kinds of rod theory with bending- 
induced softening. First, let us assume that the bend- 
ing moment depends linearly on the curvature for small 
curvature but not on the curvature for large curvature, 
which is expressed as 

M ={ ^ < Kc -* (29) 

[ klK c , [K > K c ) 

where k\ and k c are the elastic bending rigidity and the 
critical curvature, respectively. Eq. (|29]) describes the 
bending-induced softening relation of the first kind which 
is depicted in Fig. Ela) . The corresponding free energy 
density can be expressed as 

G = 7 + (fci/2)^ 2 - (« - k c ) 2 H{k - k )}, (30) 

where Tt(.) is the Heaviside step function. The above 
form has been employed by Yan et al. to inves- 
tigate the loop formation mechanism and probability 
of short DNA rings ^ We conjecture that this model 
could solve the paradox in the experiment on the ring 



Consider a rod divided into two parts at s = L c : one 
part (s < L c ) has curvatures less than k c another one 
larger than k c . In terms of the variational method in 
Appendix [XI we can derive the equations describing the 
rod as 

fci(2«;" - 2kt 2 + k 3 ) - 2 7 k = {s<L c ), (31) 
kt' + 2k't = (s < L c ), (32) 
k\K c (K c K — 2t 2 ) — 27K = (s > L c ), (33) 
r'=0 (s>L c ). (34) 

At the divided point s = L c , we have the joint conditions 
as 

K_ = K + = K c , (35) 

k'_ = 0, (36) 
t-=t+, (37) 

where (.)_ and (.)+ represent the values of (.) at the left 
and right sides of s = L c . 

Secondly, let us assume that the bending moment de- 
pends linearly on the curvature for small curvature but 
weaker linearly on the curvature for large curvature, 
which is expressed as 

M = I ( K ^ Kc ) (QQ) 

\ k 2 (K - k c ) + kiK c , (k>k c ) 

where k\ > fc 2 are the elastic bending rigidities while k c 
is the critical curvature. Eq. ([38]) describes the bending- 
induced softening relation of the second kind which is 
depicted in Fig. E^b) . The corresponding free energy den- 
sity can be expressed as 

G = 1 +(k 1 /2)K 2 +[(k 2 -k 1 )/2}(K-K c ) 2 H(K~K ), (39) 
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Consider a rod divided into two parts at s = L c : one 
part (s < L c ) has curvatures less than n c another one 
larger than k c . In terms of the variational method in 
Appendix [XJ we can derive the equations describing the 
rod as 

/ci(2k"-2kt 2 + k 3 )-27k = (s < L c ), (40) 
kt'+2k't = Q (s<L c ), (41) 
2k 2 K N + [k 2 (K - k c ) + kiK c ](K 2 ~ 2t 2 ) - 27K 

+(k 2 - ki){n - k c )kk c = (s>L c ), (42) 
k 2 {T'n + 2k! t) + (fci - k 2 )n c T' = (s > L c ). (43) 

At the divided point s = L Cl we have the joint condi- 
tions as 



/ci(k_ - K c ) = k 2 (K + - K c ) 
k\K.'_ = k 2 K' + 



T_ = T. 
.2 



/ci(k_ - k + ) = (/c 2 - _ K c) z 



(44) 
(45) 
(46) 
(47) 



Obviously, the above equations (f40 |) -(|47 |) degenerate 
into Eqs. dSID dSZJ) if fc 2 = and into Eqs. flll])-([ll]) if 
k 2 = ki . 



C. Elastic theory of 2D continua 

A 2D continuum can be simplified as a surface as 
shown in Fig. [|5J At each point, we can select a frame 
{ei, e 2 , e^}. A pressure p is loaded on the surface in the 
inverse direction of the normal vector 63. Let us cut a 
region enclosed in any curve C from the surface, t is the 
tangent vector at point of curve C. b is normal to t and 
in the tangent plane. The force and moment per length 
performed by the other region on curve C are denoted as 
f and m, respectively. Through Newton's laws, the force 
and moment balance conditions are obtained as 



f ds — I pe$ dA = 0, 



(48) 



® mds + (p rxfds— / rx pe 3 dA = 0, (49) 
Jc Jc J 

where ds and dA are the arc length element of curve 
C and area element of the region enclosed in curve C, 
respectively. 

Define two second order tensors 6 and DJl such that 



e-b = f, OTb = m. 



(50) 



These two tensors can be called as stress tensor and bend- 
ing moment tensor, respectively. Using the Stokes' theo- 
rem, we can derive 



/ 



(div 6 - pe 3 ) dA = 0, (51) 
(div0Jl + ei x ©i +e 2 x & 2 )dA = 0. (52) 




FIG. 



Force and moment in a 2D continuum. 



where Si = <S ■ B\ and & 2 = 6 ■ e 2 . Since the integral 
is performed on the region enclosed in an arbitrary curve 
C, from the above two equations we obtain the force and 
moment balance conditions of 2D continua as: 

div & = pe 3 , (53) 
div aji = 6 1 xei+6 2 xe 2 . (54) 

The above two equations are equivalent to Eq. (25) in 
Ref.[H, and Eqs. (28) and (57) in Rei.M- Eqs. ^ and 
(|54p with some complement constitutive relations form 
the fundamental equations of 2D continua. 



1. Fluid membranes 

A fluid membrane is a 2D isotropic continuum which 
cannot withstand in-plane shear strain. Generally, we 
assume that the fluid is incompressible. The free energy 
density, G, of fluid membranes should be invariant un- 
der the in-plane coordinate transformation. In terms of 
the surface theory, there are only two fundamental geo- 
metric invariants: the mean curvature 2H and gaussian 
curvature K. Thus the free energy density should be a 
function of 2H and K, that is, 



G = G(2H,K). 



(55) 



The free energy of a closed fluid membrane can be 
expressed as 



T = \ GdA+p / dV, 



(56) 



where dA is the area element of the membrane and dV 
is the volume element enclosed in the membrane, p is 
the osmotic pressure, the pressure difference between the 
outer and inner side of the membrane. The general Euler- 
Lagrange equation of free energy (|56p can be derived 
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through the variational method shown in Appendix [B] 
as 

p - 2HG + (V 2 /2 + 2H 2 - K)(dG/dH) 

+ (V -V + 2KH)(dG/dK) = Q. (57) 

As we known, the above equation has been derived by 
several authors such as Ou-Yang et all 49 ' 58 and Giaquinta 
et al~ coming from different research fields. It is re- 
cently employed to investigate the modified Korteweg-de 
Vries surfaces^ Here V • V can be called as the Laplace 
operator of the second class which is also fully discussed 
by Zhang and Xu£i 

We emphasize that ([57]) can be also derived from the 
bottom-up method, Eqs. (|53[) and (|54[) combining a com- 
plement constitutive relation 

SOI = (G 6 /2)(eiei - e 2 e 2 ) - G a e 2ei + G c e ie2 , (58) 

where G a , Gb, and G c represent the partial derivatives of 
G with respect to a, b, and c, respectively. Here a, &, and 
c are the components of the curvature tensor in Eq. {5]). 
To illuminate this point, we consider an example in which 
the free energy density is taken as G — k c (2H) 2 +X, where 
k c and A are the bending modulus and surface tension of 
the fluid membrane. It follows that Tl = 2k c H(eie2 — 
e 2 ei) from Eq. ([55)) . Substituting it into Eqs. ((53]) and 
([S"4"]) . we can derive 

p - 2XH + 4k c H(H 2 -K) + 2k c V 2 H = 0, (59) 

which is the same as the result obtained directly from 
(|57p . Simultaneously, we have the stress components 

6i = (2H 2 - 2aH + X)e 1 - 2bHe 2 - 2Hie 3 , (60) 
e 2 = -2bHe 1 + (2H 2 - 2cH + A)e 2 - 2H 2 e 3 , (61) 

where Hi and H 2 are the directional derivatives of H 
respect to ei and e 2 . These equations have been also de- 
rived by Capovilla and Guven^i from which we seem to 
arrive at a paradox for fluid membranes: we have men- 
tioned that fluid membranes cannot withstand in-plane 
shear strain, however Eqs. (|6"0|) and (|6"Tj) reveals shear 
stress still exhibits in non-spherical vesicles. 

2. Solid shells 

A solid shell is a 2D isotropic continuum which can 
endure both bending and in-plane shear strain. The free 
energy density, G, of solid shells should be invariant un- 
der the in-plane coordinate transformation. There are 
only two fundamental geometric invariants, 2H and K , 
and two fundamental strain invariants: the trace, 2 J, and 
the determinate, Q, of the in-plane strain tensor. Thus 
free energy density should be a function of 2H, K, 2J, 
and Q. That is, G = G(2H, K; 2 J, Q). 

If the solid shell has no initial strains and consists of 
materials distributing symmetrically with regard to the 



middle surface of the shell, we can expand G up to the 
second order terms of curvatures and strains as 

G = (k c /2)(2H) 2 -kK+ {k d /2){23) 2 - kQ, (62) 

where k c and k are the bending moduli while kd and k 
are the in-plane rigidity moduli. The theory based on the 
above free energy density is called Kirchhoff 's linear shell 
theory. 2 Especially, if the shell consists of 3D isotropic 
materials, we have 

k c = Y7i 3 /12(l - v 2 ), (63) 
k d = Yh/{l-u 2 ), (64) 
k/k c = k/k d = {l-u), (65) 

where Y and v are the Young's modulus and Poisson 
ratio while h is the thickness of the shells 

For a closed shell, its free energy is expressed as 
Eq. (|56"|) with G in Eq. (f6"2"| . Of course, we can obtain 
the equations of in-plane strains and shapes through the 
variational method in Appendix B. The final results are 
the same as those obtained from Eqs. and (|S"4")) with 
a complement constitutive relations (|58p and 

6 = 6 i + 6 f (66) 

with 

& = G £ll e 1 ei + (G ei2 /2)(e 1 e 2 + e 2 e 1 ) + G e22 e 2 e 2 , (67) 

where G £ll , G £l2 , and G £22 represent the partial deriva- 
tives of G with respect to en, s\2, and e 22 , the com- 
ponents of the in-plane strain tensor €. Substituting 
Eq. {55} into Eqs. {55]) and (|67|). and then employing 
Eqs. {53} and {52]), we obtain 

d(6 l n w 2 - 6 l 12 ^i) - (S 21 w 2 - &\ 2 u) X ) A w 21 = 0, (68) 
d(& 21 u 2 - & 22 uji) - {6 l n Lo 2 - ©i 2 ^i) A cj 12 = 0, (69) 

and 

p+2k c [2H{H 2 -K)+V 2 H]~<i{k d ~~k)JH-~m : £ = 0, (70) 

where 6' u = (2k d J ~ ke 22 ), &\ 2 = & 21 — kei 2 , and 
& 22 = (2kdJ — ken) are the components of tensor (5\ 91 
is the curvature tensor related to Eq. {5]). The above 
equations (|68 ]) - (|70]) describe the in-plane strains and 
shapes of solid shells at equilibrium state. The sim- 
ilar equations and the corresponding dynamics forms 
have been derived through the variational method in 
Refs. |4^ and respectively, with the aid of moving 
frame method. 

The above equations ([68]) and ([69]) can be written as 
one vector equation by introducing a displacement vector 
u = wiei+M 2 e 2 -|-U3e3, which is related to two invariants 
23 and Q of the in-plane strain tensor as 

23 = divu-2iftt 3 (71) 
2Q = (divu- 2Hu 3 ) 2 + (l/2)(curlu) 2 - (^u) 2 , (72) 
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where <)u = Vu — e3(e3 ■ Vu) is the in-plane part of Vu. 
Using the new variable u, Eqs. (p5|) and (fBU)) can be be 
written as 

(k-2k d )V{divu-2Hu 3 )-k(<> 2 u+Ku+Vu 3 ) = 0, (73) 

where u and ^ 2 u are the in-plane components of u and 
div(<)u), respectively. V is called the gradient opera- 
tor of the second class, which is shown in our previous 
work42 In particular, H, K, Vi<3 vanish and (} 2 degen- 
erates into V 2 for a flat manifold. Then the above equa- 
tion degenerates into the Cauchy equation^ in 2D plane. 
Thus Eq. ([73]) can be regarded as the Cauchy equation 
in a curved surface. 



III. APPLICATION OF ELASTIC THEORY IN 
BIO-STRUCTURES 

In the above section, we have described fundamentals 
of geometric and elastic theory on low-dimensional con- 
tinua. Can this theory be applied to the bio-structures, 
such as DNA and cell membranes, and so on? DNA is a 
long chain macromolecule which may be described as an 
elastic rod. A cell membrane is a thin structure whose 
thickness and the size of the microscopic components are 
so much smaller than its lateral dimension that it can 
be regarded as a 2D continuum phenomenologically. We 
will discuss the application of the above elastic theory in 
short DNA rings, lipid membranes and cell membranes 
in this section. 



A. Short DNA ring 

DNA is a double helical structure whose diameter is 
about 2.5 nm. Its bending rigidity, described as the per- 
sistence length l p , is about 50 nm (150 bp) at the room 
temperature. The normal DNA is usually flexible enough 
because its length is so much larger than l p that the fluc- 
tuations are quite evident. Thus the rod theory cannot 
directly be applied to the normal DNA. The statistical 
theory combining the rod theory is required ) 63 ! which 
is out of our topic in this review. However, there is a 
special kind of short DNA ring a 65 i 66 ' 67 which are in the 
length scale of l p so that the fluctuation effect can be 
neglected. The diameter is still much smaller than the 
total length. Thus the rod theory mentioned in Sec. Ill Bl 
is expected to be available for this kind of DNA rings. 

Han et al. have used AFM to observe DNA rings con- 
sisting of several segments connected by kinks in the pres- 
ence of Zn 2+ ionsi 65 ' 66 Zhao et al. have analyzed the 
mechanism of this kink instability based on Hclfrich rod 
theory^ Their main ideas are sketched as follows. First, 
a circle is a solution to Eqs. (|2T|) and (f2"5|) . Next, through 
analyzing the stability of the cycle, it is found that, for 
the given elastic constants, there exists a critical radius 
above which DNA circles will be instable. This prediction 
is in good agreement with the experiments^ 6 ^ 6 - where 



kink deformations were observed in DNA rings of 168 
bp but not 126 bp. Above some thresholds of the chi- 
ral modulus, k 3 in Eq. (|2l)|). the DNA circles turn into 
elliptical, triangular, square, or other polygonal shapes, 
respectively. This fact agrees with the experiments if £3 
is positively correlated to the condensation of Zn 2+ ions. 

Interestingly, Zhou and Ou-Yang proposed another in- 
terpretation based on the dynamic instability of Kirch- 
hoff rod theory^ with R2 = K3 = in Eq. (|10p. Their 
result is the same as that obtained directly from the first 
and second order variations of the free energy. We deal 
with the latter scenario. First, r = 0, <f) = 0, and n = 1/R 
satisfy Eqs. (f23]) - ((2"5"|) derived from the first order vari- 
ation of the free energy. That is, a planar circle with 
radius R is an equilibrium configuration. Next, through 
the second order variation of the free energy, we can ob- 
tain the characteristic function describing the stability of 
the circle 

9c {R) = r\ - (1 -T^/R-Tn 2 /R 2 < (74) 

where n > 1 is an arbitrary integer and T = k$/k\. From 
the above inequality, we obtain the critical radius 

r c = st/k! [r - 1 + y/(r - 1) 2 + i6r] , (75) 

above which the circle is instable. If only the presence of 
Zn 2+ ions tunes the values of T and fc 3 /fc 2 such that R c is 
in the range between 63/7T (bp) and 84 /ir (bp), the above 
result is also in agreement with the experiments ) 65 ! 66 
where kink deformations were observed in DNA rings of 
168 bp but not 126 bp. 




FIG. 7: A possible configuration of a short DNA ring. 

In Sec. IIIB| we also mention the theory of bending- 
soften rod. Can this theory also provide an interpretation 
to the experiments? Let us consider the bending-soften 
rod theory of the first kind whose free energy density is 
expressed as Eq. (|30]) . When the radius R of the ring is 
smaller than 1/k c , any small perturbation will increase 
the free energy. If R > l//c c , the ring might transform 
into the fictitious configuration shown in Fig. [7] which 
consists of four arcs AB, BC, CD, DA with the radius R\ 
and i? 2 - To see conveniently, the joint points are marked 
as small cycles in the figure. Obviously, R2 < R < Ri- 
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Through simple calculations, we find that the fictitious 
configuration is energetically less favorable than the per- 
fect ring with radius R. Therefore, this coarse analysis 
reveals that the theory of bending-soften rod cannot ex- 
plain the experiments. 



B. Lipid membrane 

Lipids are dominant composition of cell membranes. 
Most of lipid molecules have a polar hydrophilic head 
group and two hydrophobic hydrocarbon tails. When a 
quantity of lipid molecules disperse in water, they will 
assemble themselves into a bilayer vesicle as depicted 
in Fig. [HI in which the hydrophilic heads shield the hy- 
drophobic tails from the water surroundings because of 
the hydrophobic forces. This self-assembly process has 
been numerically investigated by Lipowsky et aZ . 68 ' 69 i 70 
and Noguchi et ali 7 ^ through molecular dynamics simu- 
lation based on coarse-grained model or meshless mem- 
brane model. 
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FIG. 8: A lipid bilayer vesicle. 

The thickness of the lipid bilayer and the size of single 
lipid molecules are much smaller than the scale of the 
whole lipid bilayer. Additionally, at the physiological 
temperature, the lipid bilayer is usually at the nematic 
state where the hydrocarbon chains of the lipid molecules 
are roughly perpendicular to the bilayer surface. Thus 
the bilayer can be regarded as a 2D fluid membrane whose 
free energy density is expressed as Eq. (|55|) . Expanding 
it up to the second order terms of curvatures, we obtain 
the Helfrich's form)22 



G H = (k c /2)(2H + c ) 2 - kK + A, 



(76) 



where k c and k are the bending moduli of the lipid bi- 
layer. We emphasize that the minus sign before k in 
Eq. (|T6"|) is opposite to Helfrich's convention. A is the 
surface tension of the bilayer. cq is called the sponta- 
neous curvature that reflects asymmetric factors between 
two sides of the bilayer, including the lipid distribution, 
the chemical environment, and so on. k c is about 20 T 
for lipid bilayers, where the Boltzmann factor is set to 1 



and T the room temperature, from which the persistence 
length of lipid bilayers is estimated about 10 /xmj 32 i 33 
In this section we only consider the size of lipid bilayers 
smaller than 10 /im so that the fluctuation effect on the 
shape of lipid bilayers can be neglected. The model based 
on Eq. f76")) is called spontaneous curvature model. We 
still remind gentle readers to note the two similar non- 
local models — the bilayer-coupling mode l 72 1 73 and the 
area difference model^ although we will not touch them 
in the present review. 



1. Closed vesicles 

The free energy of a lipid vesicle under the osmotic 
pressure p (the outer pressure minus the inner one) can 
be written as Eq. (p)6"|) with G = Gh being Helfrich's 
form (|76f . Substituting ([76)) into Eq. ([57)) , we can obtain 
the shape equation of lipid vesicles; 75 ' 76 

p-2\H + k c (2H + co)(2H 2 -c H-2K) + 2k c V 2 H = 0. (77) 

This equation is the fourth order nonlinear equation. It 
is not easy to find its special solutions. We have known 
three typical analytical solutions: sphere^ 5 , torus ) 77 ' 78 
and biconcave discoid shaped 

For a sphere with radius R, we have H = —1/R and 
K = 1/R 2 . Substituting them into ([77]), we arrive at 



pR 2 + 2XR - k c c a (2 - c a R) = 0. 



(78) 



This equation gives the sphere radius under the osmotic 
pressure p. 

A torus is a revolution surface generated by a circle 
with radius p rotating around an axis in the same plane 
of the circle. The revolution radius r should be larger 
than p. A point in the torus can be expressed as a vector 
{ (r + p cos ip) cos 9, (r + p cos <p) sin 9, p sin ip} . Through 
simple calculations, we have 2H = —(r + 2p cos ip)/p(r + 
pcosip), K = cos f/p(r + pcosip). Substituting them 
into Eq. l[77jl. we derive 



4Xp 2 + 2Pp 3 )/g 3 } cos 3 ip 
lOXp 2 + 6Pp 3 )/g 2 ]cos 2 ip 
8Xp 2 + 6Pp 3 )/g]cosip 



[(2k c clp 2 - 4k c c p 
+ [(5fc c CoP 2 - 8k c c p 
+ [(4k c clp 2 - ik c c p 
+2k c /g 2 + k c (c 2 p 2 - 



I) 



2(Pp + X)p 2 



0. 



(79) 



where g — r/p. If g is finite, then Eq. f79|) holds if and 
only if the coefficients of {1, cos ip, cos 2 <p, cos 3 ip} vanish. 
It follows 2Xp 2 = fc c co/o(4 — Cop), Pp 3 = —2k c pco and 
g = \/2. 77 That is, there exists a lipid torus with the ratio 
of its two generated radii being \f2, which was confirmed 
in the experiment [80l 

To describe the solution of biconcave discoid shape, we 
write the shape equation (|77|) under the axisymmetric 
condition. If a planar curve z = z(p) revolves around the 
z-axis, an axisymmetric surface is formed. Each point 
on the surface is expressed as r = {p cos 1^9, psin<p, z(p)}. 
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Denote ip = arctan((iz/dp) and ^ = sin?/;. Then Eq. (1771) 
is transformed into^ 1 - 



+ Co 



1 - ^ 2 



+ JL = o, (80) 



where the prime represents the derivative with respect to 
p\ This equation is called the shape equation of axisym- 
metric lipid vesicles. Its first integral, group structure 
and corresponding Hamilton's equations are investigated 
by Zheng and Lhi)22 Xu and Ou-Yangj^ 3 - and Capovilla 
et aL s4 i 85 respectively. 




P B p 

FIG. 9: A quarter outline of the biconcave surface. 

It is easy to verify that 'F = sin?/) = —cophi(p/ pg) 
with a constant ps is a solution to Eq. (|80p if p and A 
are vanishing. For < cqPb < e, the parameter equation 



sin ip = -c p\n(p/ p B ) 
z = zq + Jq tan ipdp 



(81) 



corresponds to a curve shown in Fig. [HI A biconcave 
discoid surface will be achieved when this curve revolves 
around z-axis and then reflects concerning the horizontal 
plane. The above equation ([8Tj) can give a good explana- 
tion to the shape of human red blood cell under normal 
physiological conditions i 7 - 9 - If CqPb is out of the range be- 
tween and e, Eq. (18 lj) corresponds to a prolate ellipsoid 
or other self-intersecting surfaces^ 

In the purely mathematical viewpoint, there are also 
the other solutions to Eq. ((77)) such as cylinder, constant 
mean curvature surface, periodic undulation surface^ 
pearling tubule,— and so on . 34 ' 89 However, It is a pity 
that they are open surfaces and do not correspond to 
truly closed vesicles. 

As mentioned above, it is fairly difficult to find the 
analytical solution to Eq. ([77)) . Thus we appreciate the 
applications of numerical methods to find the equilib- 
rium shapes of closed vesicles. Two kinds of typical 
numerical frameworks are usually employed. The first 
one is to use 'Surface Evolver', a software package de- 
veloped by Brakke^ to find the configurations minimiz- 
ing the free energy under some constraints . 91 ' 92 ! 93 ' 94 The 
second one is based on the phase field formulation of 
Helfrich's free energy density ([75)) and diffusive interface 
approximation ! 95 ' 96 ' 97 ' 98 The above numerical methods 



can obtain lipid vesicles with different shapes either ax- 
isymmetric or asymmetric. Additionally, the finite ele- 
ment method might be a potential method although very 
sparse literature— treats lipid bilayers by using it. 



2. Stability of closed vesicles 

When the osmotic pressure is beyond some threshold, 
a closed vesicle will lose its stability and change its shape 
abruptly. The threshold is called the critical pressure. To 
obtain it, one should calculate the second order variation 
of the free energy (J56j with G being Helfrich's form (|76|) , 
which has been dealt with in the general case as ; 49 ' 100 



5 2 T = J fc c [(V 2 fi 3 ) 2 + (2H + c )V(2#fi 3 ) ■ VQ 3 ]dA 

+ J [4k c (2H 2 - K) 2 + k c K{cl - AH 2 ) + 2XK - 2Hp\Q. 2 3 dA 

+ J [k c (UH 2 + 2c H -4K - c Q /2) - X]a a V 2 n 3 dA 

J {2H + c )[Vfi 3 • Vfi 3 + 2^3 V ■ VQ 3 ]dA, (82) 



-2k 



where ^3 is an arbitrary small out-of-plane displacement 
and the operator V is the gradient operator of the second 
class^ 9 - 

Here we will mention two results for special configura- 
tions. 

First, let us consider a lipid sphere that satisfies 
Eq. (|78[) . On the sphere, the function can be ex- 
panded by the spherical harmonic functions Y) m as = 



1=0 l~mi=-l a lr, 



Yi m . Substituting it into Eq. 



derive 8 2 T = (R/2) £ i m \a lm \ 2 [1(1 + 1) - 2}{2k c [l(l + 
- — p}, From which we can obtain the critical 



l)-c R]/R 3 
pressure^ 



Pc = 2fc c (6 - c R)/R 3 



(83) 



If p < Pc, 5 2 T > for any |o/ m |; on the contrary, 8 2 T 
can be negative for the special selection of |aj m |. The 
above equation depends also on cq. If Co > G/R, then p c 
is negative, which reveals that a sphere vesicle is always 
instable for large enough cq. 

Next, let us still regard a long enough lipid tubule as 
a closed vesicle. Denoted its radius as p. From Eq. ([77]) 
we have 



(k c /2)(l/p 2 -c 2 )-pp = X. 



(84) 



On the cylindrical surface, can be expanded as Fourier 
series Q3 = 53z=-oo a; exp(iZ0). Substituting it into 
Eq. ([g2"| and combining Eq. (|53)l. we derive 8 2 T = 



J2i=-oo \ a l\ v> — — — p]i From which we 

can obtain the critical pressure 

3k c /p 3 . 



Pc 



(85) 



If p < p c , 8 2 T > for any |a;|; on the contrary, 8 2 T can 
be negative for the special selection of |aj|. 
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3. Open vesicles with free edges 

The opening-up process of lipid vesicles by talin, a pro- 
tein, has recently been observe d 101 ' 102 which pushes us 
to study the equilibrium equation and boundary condi- 
tions of lipid vesicles with free exposed edges. Capovilla 
et al. have addressed this problem and given the equi- 
librium equation and boundary conditions Inspired 
by the talk "moving frame method" of Chern,-i2^ we in- 
troduce exterior differential form to deal with the varia- 
tional problem on open surface and obtain concisely the 
shape equation and boundary conditions of open lipid 
vesicles Numerical solution to the shape equation and 
boundary conditions with relaxed method can explain 
the experimental results very wellJ^? A quantity of open 
vesicles with free edges have also been obtained numer- 
ically by Wang and DuiSI with the phase field method. 
Here we will not further discussed the dynamical opening 
process of the vesicles, which has been recently investi- 
gated by Kaga and Ohta.*^ 




FIG. 10: An open surface with boundary curve C. 

We regard an open lipid vesicle with a free edge as a 
smooth surface with a boundary curve C, as shown in 
Fig. [101 t is the tangent vector of the curve C. b, in the 
tangent plane of the surface, is perpendicular to t and 
points to the opposite side that the surface located in. 
The free energy of the open lipid vesicle is written as 

T = J G H dA + j <j> ds, (86) 

where 7 represents the line tension of the edge and Gh 
has the Helfrich's form ([TH)) . 

The first order variation of T gives the shape equation 

k c (2H+c )(2H 2 -c H-2K)-2\H+2k c V 2 H = 0, (87) 

and the boundary conditions aSf^S^ 

[k c {2H + c )-kk n ] c = 0, (88) 
[2k c dH/dh + 7 fc„ - kr' g ] c = 0, (89) 
[G H + jk g ] c = 0, (90) 

where fc„ and k g are normal curvature and geodesic cur- 
vature of the boundary curve C. r g is the derivative of 
geodesic torsion r g with respect to the arc length of curve 
C . The mechanical meanings of the above four equations 



are as follows: Eq. ((87)) is the normal force balance equa- 
tion of the membrane; Eq. (|88p is the moment balance 
equation of points in curve C around the direction of t; 
Eq. ([8"9")) is the force balance equation of points in curve 
C along the normal direction of surface; and Eq. ([90]) 
is the force balance equation of points in curve C along 
the direction of b. It is necessary to emphasize that the 
boundary conditions are available for open vesicles with 
more than one free edge because the edge in our deriva- 
tion is a gener al one. 

In Ref. Il05l we have shown two analytical solutions to 
above equations ([5T[) - ([9T)|) : One is a cup- like membrane 
and another is the central part of a torus. Several numeri- 
cal solutions to these equations are obtained by Umeda et 
alr^. Their results reveal that the line tension 7 induced 
by talin correlates negatively with the concentration of 
talin, which is in agreement with the experimental result 
that the hole of vesicle is enlarged with the concentration 
of talin.i^ 1 



4- Vesicles with lipid domains 

The above discussion on open lipid vesicles with free 
edges can be extended to study a vesicle of several lipid 
components. The domains usually formed so that each 
domain contains one or two kinds of lipid molecules. 
The morphology of axisymmetric vesicles with multi- 
domains has been theoretically investigated by Jiilicher 
and Lipowsky^ ^ It is found that lipid domains facilitate 
the budding of vesicles The giant vesicles with lipid 
domains have been observed in recent experiment p^ii 
There are two kinds of lipid domains which are at the 
liquid-ordered state and liquid-disordered state, respec- 
tively. It is natural to assume that different kinds of 
domains have different bending moduli and spontaneous 
curvatures. The axisymmetric vesicles in the experiment 
can be explained with Jiilicher-Lipowsky theory through 
numerically method. Baumgart et al. have demonstrated 
that the line tension, the osmotic pressure, the relative 
bending moduli, and the spontaneous curvature have 
significant effects on the morphology of a vesicle with 
two domains being at the liquid-ordered and disordered 
states, respectively*!^ 

The asym metri c vesicles are also experimentally ob- 
served in Ref. which enlightens us to investigate the 
shape equation of each domains and the boundary con- 
ditions between domains without any axisymmetric as- 
sumptions. Let us consider a vesicle with two domains 
separated by curve C sketched in Fig. [TT] The free energy 
can be expressed asi^S- 

T = J G I H dA + J Gj^dA + j Ids+pJ dV, (91) 

where and G^j have the Helfrich from ([76]) with the 
bending moduli fc 7 , k 1 , kj , k 11 , the spontaneous curva- 
tures Cg, Cq 7 , and the surface tensions A 7 , A 77 , respec- 
tively. The integrals in the first and second terms of 
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Eq. (|91[) are performed on the domain I and II shown in 
Fig. 111! respectively. 7 is the line tension of boundary 
curve C . p is the osmotic pressure of the vesicle. 




FIG. 11: A vesicle with two domains separated by curve C. 

In terms of the physical meanings of Eqs. (j87|) - ([90|) . 
we can easily write down the shape equation of domains 

p-2\ l H +kl(2H +c )(2H 2 -<%H -2A") + 2fc*V 2 H = 0. (92) 

where the superscript i = I and II represents the physical 
quantity of lipid domains I and II, respectively. Addi- 
tionally, the boundary conditions between domains are 
as follows)42ii!a 

[ki(2H + cl) - k I c I (2H + ci 1 ) - {¥ - k u )k n ] c = 0,(93) 
[2(fc c / + k^dH/db - (F + k")r' g + jk n } c = 0, (94) 
[G 1 ~G H +jk g }c = 0, (95) 

where b is perpendicular to the boundary curve C and 
points to the side of domain II. 

As we know, there is still no any numerical result on 
asymmetric vesicles with domains directly from the a bove 
equations in the previous literature. Only in Ref. 1 1071 
Wang and Du discussed the morphology of asymmetric 
vesicles with domains through the phase field model. 

In the above theory, the detailed architecture of liquid- 
ordered and disordered phases is neglected. There 
are special lipid domains at liquid-ordered phase, so 
called rafts, which are enriched in cholesterol and 
sphingolipidsjii^ Cholesterol is a kind of chiral molecules, 
which has not been included in the above theory. Re- 
cently, a concise theory of chiral lipid membranes devel- 
oped by Tu and Seifertf^ might be extended to discuss 
the raft domains. 



5. Adhesions of Vesicles 

Cell adhesion is a complex biological process which 
controls many functions of life. It can be understood as 
a first-order wetting transition 1 ^ and might be simpli- 
fied as the adhesion of lipid vesicles. As a model, Seifert 
and Lipowsky have theoretically investigated a lipid vesi- 
cle adhering to a flat rigid substrate and found that the 
vesicle undergoes a nontrivial adhesion transition from 
the free state to the bound state, which is governed 



by the competition between the bending and adhesion 
energiesjiii Ni et al. have discussed the adhering lipid 
vesicles with free edges and the adhesion between a lipid 
tubule with a rigid substrate i 118 ' 119 A big progress on this 
topic is recently made by Guven and his coworker o 120 ' 121 
who obtain the general equations to describe the con- 
tact line between the vesicle and the rigid substrate or 
another vesicle. 



(a) (b) 

Vesicle Vesicle I 




Substrate Vesicle II 



FIG. 12: Adhesions, (a) Adhesion between a lipid vesicle and 
rigid substrate with a contact line C. (b) Adhesion between 
two lipid vesicles with a contact line C. 

The adhesion between a lipid vesicle and a rigid sub- 
strate is depicted in Fig. [I^k where the contact area is 
denoted by A. The free energy of this system is expressed 

T = J G H dA + pJ dV- WA, (96) 

where p is the osmotic pressure of the vesicle and W is the 
strength of the adhesion potential between the vesicle and 
the substrate. Gh is the free energy density of Helfrich's 
form ([75)1 . For the flat rigid substrate, a characteristic 
radius and the length scale of the vesicle are defined as 
R a = yj2k c /W and R = y/A/4%, respectively. If R < 
R a , the vesicle is a little stiffer or the attraction is relative 
weak such that A approaches to zero. Thus the vesicle 
is unbound to the substrate and this state is called the 
free state. On the contrary, the vesicle is at the bound 
state. At this state, let us take t as the tangent vector 
of the contact line C, and b perpendicular to t and in 
the common tangent plane of the lipid vesicle and the 
substrate. The absolute value of the normal curvature 
al ong b f or the point on the contact line is proven to be 
y/2W/k c for an axisymmetric vesicle adhering to the flat 
substrateAii If the rigid substrate is curved, the above 
conclusion is revised as^2i 

I «b -«b I = y/2W/k c , (97) 

where and are the normal curvatures along b for 
the points outside but near the contact line, calculated 
by using the surfaces of the vesicle and the substrate, 
respectively. 

The adhesion between two lipid vesicles is depicted in 
Fig. [T!?b . The free energy of this system is expressed 
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J G' H dA + J p I dV + J GndA + J 



P n dV-WA, 



(98) 



where p 1 and p* 



jI are the osmotic pressures of the vesi- 
cles I and II, respectively. A and W are the contact area 
and adhesion strength, respectively. G ! H and G 1 ^ are 
the Helfrich's free energy density of vesicle I and II. The 
first order variation of ([55]) gives the same shape equa- 
tion of two vesicles as (|92|) and the adhesion boundary 
conditions 



(1 
(1 



= 2W/kl 



#7*c)(«" 

JI 



A\2 
b ^h) 



2W/k 



ii 



^ - K £)/dh = 0, 



(99) 
(100) 
(101) 



where and n{f are the normal curvatures along b for 
the points outside the adhesion domain but near the con- 
tact line calculated by using the surfaces of vesicles I and 
II, respectively, is the normal curvature for the points 
inside the adhesion domain but near the contact line cal- 
culated by using the common surface of vesicles I and 
II. As we know, there is still lack of numerical solutions 
to the above equations (f9"9"|) - (|101[) in the previous litera- 
ture. Only in the recent work, Ziherl and Svetina^^ have 
investigated the adhesion between two vesicles by numer- 
ically minimizing the free energy (|98p with k l c l = k\ and 
various W. 

Is the behavior of vesicle adhesion close to that of 
cell adhesion? The cell membrane can bear shear strain 
whose adhesion behavior might be much closer to the 
adhesion between a polyelectrolyte microcapsule and the 
substrate*!^ 3 - Interestingly, beyond the threshold adhe- 
sion strength W c , the contact length scale increases in 
proportion to (W — Wc) 1 / 2 , which is the same as the be- 
havior of vesicle adhesions except the coefficient before 
(W- Wc) 1 ' 2 . 



6. A different viewpoint of surface tension 

Although the lipid bilayer cannot withstand the in- 
plane shear strain, it can still endure the in-plane com- 
pression strain. The in-plane compression modulus, kb, 
of lipid bilayers is about 0.24 N/m»i2^ Considering this 
point, we may write the free energy of a closed lipid vesi- 
cle as 



T = p dV 



G B dA 



(k b /2)(2J b ) 2 dA, (102) 



where 



G B = {k c /2)(2H + c ) 2 -kK, 



(103) 



and Jb is the in-plane compression or stretch strain. We 
emphasize that the contribution of chemical potential are 
omitted when we write the above free energy. 



The first order variation of the free energy (| 1 02|) reveals 
that 2 Jb is a constant and then 

p-2(2k b J b )H + 2k c V 2 H 

+k c (2H + c ){2H 2 - c H - 2K) = 0. (104) 

Comparing the above equation with the shape equation 
([77]) of lipid vesicles, we deduce that 



A = 2kbJb- 



(105) 



In the discussion on the stability of closed lipid vesicles, 
we have seen that the surface tensor A has no effect on the 
critical pressure. The second order variation of the free 
energy (I102j) can give the same conclusion. 5 2 [p J dV + 
J GsdA] has been shown in Eq. (|82|) with vanishing A. 
The additional term is 

S 2 J(k b /2)(2Jb) 2 dA = J fc h (divv- 2HQ 3 ) 2 dA (106) 

where v = f2iei + ^2^2 + ^363 represents the infinites- 
imal displacement vector of the vesicle surface. We can 
always select the proper deformation modes such that 
divv - 2^3 = and then 8 2 $ {k b /2){2J b ) 2 dA van- 
ishs, but S 2 [p j dV + J GsdA] is not affected. That is, 
the critical pressure is determined merely by S 2 [p J dV + 
J Grid A], which is independent on the compression mod- 
ulus of lipid bilayer kb- 



C. Cell membrane 

Cell membrane consists of lipids, proteins, and a small 
quantity of carbohydrates and so on. A simple but 
widely accepted model for cell membranes is the fluid 
mosaic modeli2^ proposed by Singer and Nicolson in 
1972. In this model, the cell membrane is considered 
as a lipid bilayer where the lipid molecules can move 
freely in the membrane surface like fluid, while the pro- 
teins are embedded in the lipid bilayer. Some pro- 
teins, so called integral membrane proteins, traverse en- 
tirely in the lipid bilayer and play the role of informa- 
tion and matter communications between the interior 
of the cell and its outer environment. The others, so 
called peripheral membrane proteins, are partially em- 
bedded in the bilayer and accomplish the other biolog- 
ical functions. Beneath the lipid membrane, the mem- 
brane skeleton, a network of proteins, links with the pro- 
teins embedded in the lipid membrane. Mature mam- 
malian and human erythrocytes (i.e., red blood cells) 
are lack of a cell nucleus. Thus they provide a good 
experimental model for studying the mechanical prop- 
erties of cell membranes . 126 ' 127 ' 128 i 129 On the theoreti- 
cal side, spontaneous curvature model,— rubber mem- 
brane model ) 30 i 130 i 131 and dual network models have 
been employed to investigate the mechanical and thermal 
fluctuation properties of erythrocyte membranes. We 
will address the elasticity and stability of composite shell 
model for cell membranes in this section. 
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1. Composite shell model of cell membranes 

A cell membrane can be simplified as a composite 
shells of lipid bilayer and membrane skeleton. The 
membrane skeleton, inside of the cell membrane, is a 
network of protein filaments as shown in Fig. 1131 The 
joint points of the network are bulk proteins embedded 
in the lipid bilayer. The whole membrane skeleton seems 
to float the sea of the lipid bilayer. It can have a global 
movement along the surface of the bilayer but the move- 
ment of the joints along the normal direction is totally 
coupling with the bilayer. In the mechanical point of 
view, the lipid bilayer can endure the bending deforma- 
tion but hardly bear the in-plane shear strain. On the 
contrary, the membrane skeleton can endure the in-plane 
shear strain but hardly bear the bending deformation. 
The composite shell overcomes the shortage of the lipid 
bilayer and the membrane skeleton. It can sustain both 
bending deformation and in-plane shear strain. 



nside of cell membrane 




Outide of cell membrane 

FIG. 13: Local schematic picture of the composite shell model 
for a cell membrane. 

The contour length of protein chain between joints 
in the membrane skeleton is about 100 nm which is 
much smaller than the size (~ 10 /im) of cell mem- 
branes. The lipid bilayer is 2D homogenous. The mem- 
brane skeleton is roughly a 2D locally hexagonal lat- 
tice. As is well known, the mechanical property of a 
2D hexagonal lattice is 2D isotropic. 134 Thus the com- 
posite shell of the lipid bilayer plus the membrane skele- 
ton can still be regarded as a 2D isotropic continuum. 
Its free energy density should be invariant under the in- 
plane coordinate transformation and can be written as 
G cm — G cm (2H, K;2J,Q). We can expand G cm up to 
the second order terms of curvatures and strains as 



G c 



G B + (k b /2)(2J b ) 



G 



sk ; 



(107) 



where Gb results mainly from the bending energy of 
the lipid bilayer, which has the form as Eq. ()103|) . 
(fe&/2)(2 Jj,) 2 is the contribution of in-plane compression 
of the lipid bilayer where kb and 2 Jb are the compres- 
sion modulus and relative area compression of the lipid 
bilayer. G s k — (k d /2)(2J) 2 — kQ is the in-plane com- 
pression and shear energy density which comes from the 
entropic elasticity of the membrane skeleton. k d and k 



are the compression and shear moduli of the membrane 
skeleton, respectively. Their values are experimentally 
determined as k d — k = 4.8 /xN/mi 129 i 135 2 J and Q are 
the trace and determinant of the stain tensor of the mem- 
brane skeleton. Because there is no in-plane coupling be- 
tween the lipid bilayer and the membrane skeleton in the 
composite shell model, thus J& for the lipid bilayer and J 
for the membrane skeleton have no local correlation. In 
the above subsection, we have mentioned that the effect 
of (kb/2)(2Jb) 2 can be replaced with the surface tension 
A = 2kbJb- Considering a closed cell membrane under 
osmotic pressure p, the free energy can be written as 



F= GcmdA + p dV. 



(108) 



Similarly to Sec. Ill CI if we define a displacement vec- 
tor u satisfying Eqs. ([7Tj) and ([72]) . we can derive the 
Euler-Lagrange equations corresponding to the free en- 
ergy (fTfjg)) as 

(k - 2fc d )V(2 J) - fc(^ 2 u + Ku + Vu 3 ) = 0, (109) 
p + 2k c [(2H + c )(2H 2 - c H - 2K) + 2V 2 H] - 2XH 
+2ff(fc-/c d )(2J)-fc£K:Vu = 0, (110) 

where u and <0> 2 u are the in-plane components of u and 
div(^)u), respectively. 9\ is the curvature tensor related 
to Eq. ([5]). V is called the gradient operator of the second 
class, which is shown in our previous work4^ 

Generally speaking, it is difficult to find the analytical 
solutions to Eqs. (|109[) and (|110|) . But we can verify that 
a spherical membrane with homogenous in-plane strains 
satisfy these equations. The radius R and the homoge- 
nous in-plane strain e should obey the following relation: 

pR 2 + 2(A + 2k d e - ke)R + k c c {c Q R - 2) = 0. (Ill) 



2. Stability of cell membranes and the function of 
membrane skeleton 



When the osmotic pressure is beyond some threshold, 
a closed cell membrane will lose its stability and change 
its shape abruptly. The threshold is called the critical 
pressure. To obtain it, one should calculate the second 
order variation of the free energy (|108[) in terms of Ap- 
pendix B. The variational result is 

8 2 t = J k c [(v 2 n 3 ) 2 + (2# + co)v(2fm 3 ) • vn 3 ]dA 

+ J [4k c (2H 2 - K) 2 + k c K(d - AH 2 ) + 2\K - 2Hp\U 2 3 dA 

+ J [k c (UH 2 + 2c H —AK — cg/2) - X\tt 3 X7 2 n 3 dA 

-2k c J (2H + co)[Vfi 3 ■ Vfi 3 + 2fi 3 V ■ X7£l 3 ]dA 

-k d J [(v- V + 2fm 3 )(divv-2im 3 )]d^ 

+ {k/2) J (curlv) 2 cL4-fc J Kv 2 dA + k J Q 3 V • vdA 
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+k J 2#fi 3 (divv - 2HQ 3 )dA -k J : VvdA, (112) 

where v = Qiei + Q2&2 + ^3^3 is the infinitesimal dis- 
placement vector of the cell membrane whose in-plane 
component is denoted as v = fiiei + ^12^2- 

In terms of the Hodge decomposed theorem^i 3 ^ v can 
be expressed by two scalar functions and x as 



v • dr = d£l + *dx, 



(113) 



where * is the Hodge star i 49 ' 136 Then we have divv = 
V 2 f2 and curlv = V 2 x- For the spherical cell membrane 
satisfying Eq. (|llip . Eq. (|112[) can be divided into two 
parts: one is 



5 2 T, = (A/2) J [(V 2 X ) 2 + (2/R 2 ) X V 2 x]dA; (114) 
another is 

S 2 F 2 = J Ql[2cok c /R 3 + p/R + {Ak d - 2k)/R 2 ]dA 
+ J n 3 \7 2 n 3 [k c c /R + 2k c /R 2 + P R/2]dA 
+ I k c (\7 2 n 3 ) 2 dA+[(4:k d -2k)/R} J n 3 \7 2 ildA 
-/.-,, I {V 2 fl) 2 dA+ (k/R 2 ) [ flV 2 fldA. (115) 



It is easy to verify that 8 2 T\ is always positive on a 
spherical surface. Then the stability of the spherical cell 
membrane is merely determined by b 2 T2- By analogy 
with our previous workji 3 ^ we can prove that 8 2 J-2 is 
also positive if 

p < p, S '* (2 ** " % + 2 -^[l(l + 1) - coR], (116) 
[k d l(l + 1) - k]R & 

for any integer I > 2. Thus the critical pressure is 

Pc = mm{ Pl (/ = 2, 3, 4, •••)}■ (117) 

Obviously, if k — 0, i.e., the effect of membrane skeleton 
vanishes in the cell membrane, p c degenerates into the 
critical pressure (183]) of a spherical lipid vesicle. 

When kk d (2k d - k)R 2 /k c (6k d - k) 2 > 1, the critical 
pressure is derived from Eqs. (|116[) and (|117[) as 



Pc = {A/R 2 )yJ{k/k d )(2k d -k)k c 



(118) 



As an example, let us consider a cell membrane with 
typical values of k = k d = 4.8/xN/m,i29 k c = 10^ 19 J, 
and R m 10 /im. Through a simply manipulation, we find 
that kk d {2k d - k)R 2 /k c {6k d - k) 2 > 1, and so Eq. pTB)) 
holds, from which we obtain the critical pressure p c — 
0.03 Pa. However, if the membrane skeleton vanishes, 
k = 0, we calculate p c = 0.001 Pa from Eqs. (|116p and 
(I117p . This example reveals a mechanical function of 



membrane skeleton: it highly enhances the stability of 
cell membranes. 

As a byproduct, Eq. (|118j) also gives the critical pres- 
sure 



Pi 



= ^4/3(1 - v*) Y{h/Rf 



(119) 



for a spherical thin solid shell of 3D isotropic materials if 
we take k c , k d , and k as Eqs. (|6"3"]) - ([63)1 . This formula is 
the same as the classic strict result obtained by Pogorelov 
from the other mcthodJ^i 



IV. APPLICATION OF ELASTIC THEORY IN 
NANO-STRUCTURES 

In the last section, we have expatiated on the applica- 
tion of Elastic theory in bio-structures. In this section, we 
will discuss whether and to what extent this theory can 
be applied to nano-structures, especially the graphitic 
structures, such as graphene and carbon nanotubes. 



(a) 




(b) 




FIG. 14: (a) Graphene. (b) Single-walled carbon nanotube. 



A. Graphene 

Graphene is a single layer of carbon atoms with a 
2D honeycomb lattice as shown in Fig. [T4"r . It has 
been a rapidly rising star in the material science and 
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condensed-matter physical since it was successfully 
cleaved from buck graphite J^- It is found that the free- 
standing graphene might be a strictly 2D atomic crystal 
which is stable under ambient conditions*^ 1 - However, 
Mermin has theoretically proved that the 2D crystalline 
order could not exist at finite temperature.- 1 ^ 2 - There are 
two possible ways to solve this paradox: (i) The graphene 
might not be a perfect 2D crystal. Recently, Meyer 
et al. have investigated the elaborate structure of sus- 
pended graphene sheets and found that the graphene 
sheets are not genuine flat*^ They also argue that the 
graphene sheets could be stabilized by the out-of-plane 
deformation in the third dimension resulting from the 
the thermal fluctuations. 143 Fasolino et al. have also ad- 
dressed the height fluctuations by means of Monte Carlo 
simulations.™ Their result at room temperature is in 
good agreement with the experiment mentioned above, 
(ii) Mermin theorem is valid for power-law potentials of 
the Lennard- Jones type while the interaction between 
nearest neighbor atoms (covalent bond) in the graphene 
might not be of this type*!^ 

To fully understand the experimental result and pos- 
sible stable mechanism in theory, we will address the 
Lenosky lattice model^ and its revised form as follows. 

1. Revised Lenosky lattice model and its continuum limit 

We start from the concise formula proposed by 
Lenosky et al. in 1992 to describe the deformation energy 
of a single layer of curved graphite 1 ^ 

(ij) 1 (j) 

+£ 2 _ Ui ' + 63 £( ni ' U ij)( n i ' U J"0- ( 12 °) 

(ij) (ij) 

The first two terms are the contributions of bond length 
and bond angle changes to the energy. The last two terms 
are the contributions from the 7r-electron resonance. In 
the first term, is the initial bond length of planar 
graphite, and is the bond length between atoms i and 
j after the deformations. In the remaining terms, Ujj is 
a unit vector pointing from atom i to its neighbor j, and 
rii is the unit vector normal to the plane determined by 
the three neighbors of atom i. The summation Ylu\ is 
taken over the three nearest neighbor atoms j to atom 
i, and taken over all the nearest neighbor atoms. 

The parameters (ei,e2, £3) = (0.96,1.29,0.05) eV were 
determined by Lenosky et al*^ through local density 
approximation. The value of eo was given by Zhou et al. 
as eo = 57 eV/A 2 through the force-constant method* 1 ^ 
In the above energy form, the second term requires 
that the energy cost due to in-plane bond angle changes is 
the same as that due to out-of-plane bond angle changes. 
However, the experiment by inelastic neutron scattering 
techniques reveals that the energy costs due to in-plane 
and out-of plane bond angle changes are quite different 



from each other To describe this effect, we revise the 
Lenosky lattice model as 

E 9 = I £>« ~ r Q f +e u J2 £ « ' + 1/2) 2 

(ij) i (j<k) 

+e ln £(£ u» ) 2 + e 2 £(1 - n, • n,), (121) 

i (j) (ij) 

where u^- = u i3 - (n, • Uy-Jn, and = n, • Uy. 
If the three nearest neighbor atoms to atom i are la- 
beled as 1,2,3, the summation 53(j<fc) is understood as 
Ei< 3 <fc<3- The second and third terms of Eq. (|12ip rep- 
resent the energy costs due to in-plane and out-of-plane 
bond angle changes, respectively. We have omitted the 
term £3j2(ij)( n i ' u ij)( n j ' u ji) relative to the original 
Lenosky model (|120p . because its contribution is very 
small in terms of the results by Lenosky et al. 

The parameters in Eq. (j!21[) are determined by fitting 
the total energy of variously perturbed configurations of 
\fl x v7 unit cell of graphite (14 atoms). The total energy 
is obtained through the first-principles calculations (the 
ABINIT package 1 ^). The calculations are carried by 
taking Troullier-Martins pseudopotentials, 150 plane- wave 
energy cutoff of 50 Hartree, and 4x4x1 Monkhorst-Pack 
k-points 151 in Brillouin-zone. The exchange-correlation 
energy are treated within the local-density approxima- 
tion in the Ceperley- Alder formic with the Perdew- 
Wang parametrizationJ^ 3 - Our result is ro — 1.41 A, 
e = 46.34 eV/A 2 , e u = 4.48 eV, e ln = 1.04 eV, and 
€2 = 1.24 eV. The value of eo is a little smaller than that 
obtained by Zhou et al*^ from force constant method. 
The values of ei„ and e-i are very close to those of ei 
and €2 obtained by Lenosky et alJ^- from local density 
approximation. The key reason is that the main energy 
contribution in the configurations discussed by Lenosky 
et al. comes from the third and fourth term in Eq. (|12ip . 

Now let us derive the continuum limit form of the re- 
vised Lenosky lattice model (| 12 1[) by analogy with the 
method in our previous work i 42 i 44 Now consider a curved 
graphene and take a fictitious smooth surface such that 
all carbon atoms are on that surface. The in-plane 

stain can be expressed as (£,; = ( £l1 £l2 ) in the lo- 

\ £12 £22 J 

cal frame {ei, 62,63} at atom i. The bond vector 
from atom i to its neighbor j after the deformations and 
the initial bond vector r°- before the deformations satisfy 
Yij = (1+ (Si) • r?-, where I is the unit matrix. The initial 
bond vectors can be expanded to the order of 0(r^K 2 ) 
asil 

4 = (! - *o^/6)r T J - + (« J r J r3/6)B J 

+[r Kj/2 + (rl/6)dK j /ds]r 'Nj, (122) 

where j=l, 2, 3 denote three sp 2 -bond curves from atom 
i to one of its three neighbor atoms j on the graphene sur- 
face. The symbols Tj , Nj , and represent the unit tan- 
gential, normal, and binormal vectors of the bond curve 
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from i-atom to j-atom, which satisfy the Frenent theorem 
Eq. (fl]). k, t refer to the curvature and torsion while s is 
the arc-length parameter along the bond curve. Assume 
the sp 2 -bond along the geodesic curve of the graphcnc 
surface. The vectors Tj and Bj can be expressed by 
Tj = cos^ei + sm6j&2 and Bj = — sin^ei + cos6>je2, 
where 9j is the rotating angle from ei to Tj . We have the 
expressions of uy = ?ij/rij and = Nj with r»j = |r.y| 
for the deformed graphene. Then Eq. (|121|) is trans- 
formed into the continuum limit up to the second-order 
magnitudes of en, £221 £12 and r n as 



transformed into 



E g = 



h { 2H) 2 ~kK+^-(2J) 2 ~~kQ 



with four parameters 

k c = (9ei„ + 6e 2 )r 2 /8fi 0j 
k = 3e 2 r 2 1 /4fi , 
k d = 9{e rl+3e lt )/16fl (h 
~k = 3{e rt+9e lt )/8n , 



dA, (123) 



(124) 
(125) 
(126) 
(127) 



where i7o = 3-v/3Vo/4 is the occupied area per atom. 
The continuum form (|123p has first derived in our previ- 
ous work which is, in fact, the natural conclusion of 
the symmetry of graphene»i^ The curved graphene com- 
prises a lot of hexagons which has approximately local 
hexagonal symmetry. In fact, 2D structures with hexag- 
onal symmetry are 2D isotropic.- 134 Thus the elasticity 
of the graphcnc can be reasonably described by the shell 
theory of 2D isotropic materials mentioned in Sec. Ill CI 
and so its energy has the form of Eq. (|123p . We also 
notice that a flaw in the coefficient before e\ in the ex- 
pression of k c in our previous work [3 

Using the values of tq, eit, ei n , and e 2 obtained from 
the first-principles calculations, we have k c = 1.62 eV, 
k = 0.72 eV, k d = 22.97 eV/A 2 , and k = 19.19 eV/A 2 . 
Because the results of first-principles calculation are 
applicable for zero temperature, only the results de- 
rived from the experiments at low temperature can be 
used as reference values to compared with them. The 
value of k c is close to the value 1.77 eV estimated by 
Komats u 155 i 156 i 157 at low temperature (less than 60 K). 
The value k/kg — 0.83 is quite close to the experimen- 
tal value 0.8 derived from the in-plane elastic constants of 
graphite* 1 ^ The elastic properties of graphene can be de- 
scribed by Eq. (|123p with four parameters k c , k, kd, and 
k, where the energy density is the same as Eq. (j6"2"|) . the 
free energy density of solid shell with 2D isotropic materi- 
als. Since k/k c = 0.44 is much smaller than k/kd = 0.83, 
the graphene cannot be re gard ed as a solid shell with 3D 
isotropic materials as Ref. |44|. 



E g = Qtc/2) / (V'wY d 2 x, 



(128) 



where x = (x\, X2) represents the point on the graphcnc 
plane before deformations. 
Adopting the Fourier series 



w ( x ) = ^2 S exp{iq • x}, 



(129) 



with q = (21tt/L, 2nir/L), we transform Eq. (|128p into 
£ s = (fc c /2)^qV q | 2 , (130) 



and then the corresponding partition function is derived 
as 

Z = Jlldw^ eM-E g /T) = Y[ V2vrT/fc c q 4 , (131) 



where the Boltzmann constant has been set to 1. It fol- 
lows that the equipartition theorem: 



4|„r, |2\ 



((fc c /2)q 4 |w 



-TdlnZ/dlnq 4 = T/2, (132) 



where (.) represents the ensemble average. The above 
equation is equivalent to 



(K| 2 ) = T/fc c q 4 - 
Similarly, (w 2 ) is derived as 



E 



\\ 2 ) 



TL 2 



L 2 16tt 4 /c c ^ (I 2 + n 



y L 



(133) 



(134) 



Through simply numerical manipulations, we have 1 ^ 
, „ TL 2 



{w 



150k ' 



(135) 



for the graphene co ntain s more than 100 atoms. 

In terms of Ref. Il57l . we estimate k c w 0.46 eV at 
T = 300 K. Substituting it into Eq. (fT33|) and taking L = 



25 nm as the experiment Il43l . we have y/ (w 2 ) w 0.5 nm. 
This value is a little smaller than the largest out-of-plane 
deformation 1 nm in the experiment. However, they are 
consistent with each other because 0.5 nm is the mean 
square value which should be smaller than the largest 
out-of-plane deformation in the experiment. 



2. Intrinsic roughening in graphene at temperature T 

Let us consider the freely suspended graphene which 
is almost a flat layer with the area L 2 . The small out-of- 
plane displacement is denoted by w. The energy (| 1 23[) is 



B. Carbon nanotube 

There are two kinds of carbon nanotubes: single- and 
multi- walled carbon nanotubes, which are synthesized in 
the last decade of 20 century 160 ! 161 Simply speaking, a 
single-walled carbon nanotube (SWNT) can be regarded 
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as a seamless cylinder wrapped up from a graphitic sheet, 
as shown in Fig. 114b . whose diameter is in nanometer 
scale and length from tens of nanometers to several mi- 
crometers if we ignore its two end caps. A multi-walled 
carbon nanotube (MWNT) consists of a series of coaxial 
SWNTs with layer distance about 3.4 A. 

SWNTs can be expressed as a pair of integers (n,m), 
so called index, in terms of the wrapping rule. They 
are divided into two classes: achiral tubes if n = m or 
nm = and chiral tubes for others.™ The electronic 
properties of SWNTs depend sensitively on the index)^ 
they are metallic if n — m is multiple of 3, else semi- 
conductor. SWNTs also possess many novel mechani- 
cal properties j 1 ^ 4 - in particular high stiffness and axial 
strength, which are not sensitive to the tube diameters 
and chirality. MWNTs have the similar mechanical prop- 
erties to SWNTs . 165 ' 166 In this section, we will review the 
theoretical and numerical results on the elastic proper- 
ties of carbon nanotubes, and then discuss how the low- 
dimensional elastic theory mentioned in Sec. [TT] can be 
applied in carbon nanotubes. 



1. General review on the elasticity of carbon nanotubes 



TABLE I: Young's modulus Y (unit in TPa), Poisson ratio v 
and effective thickness h (unit in A). (MD = mofecufar dy- 
namics; TB = tight-binding; SM = structure mechanics; FEM 
= finite eiement method; LDA = local density approach) 



Authors 


Y 


v 


h 


Method 


Refs. 


Yakobson et al. 


5.5 


0.19 


0.66 


MD 


40 


Tu & Ou-Yang 


4.7 


0.34 


0.75 


LDA 


44 


Kudin et al. 


3.9 


0.15 


0.89 


ab initio 


167 


Zhou et al. 


5.1 


0.24 


0.74 


TB 


168 


Vodenitcharova et al. 


4.9 




0.62 


ring theory 


169 


Pantano et al. 


4.8 


0.19 


0.75 


SM & FEM 


170.171 


Chen and Cao 


6.8 




0.80 


SM 


172 


Wang et al. 


5.1 


0.16 


0.67 


ab initio 


173 


Sears & Batra 


2.5 


0.21 


1.34 


MD 


174 


Tserpes et al. 


2.4 




1.47 


FEM 


175 


Lu 


1.0 


0.28 


3.4 


MD 


41 


Hernandez et al. 


1.2 


0.18 


3.4 


TB 


176 


Shen & Li 


1.1 


0.16 


3.4 


force- field 


177 


Li & Chou 


1.0 




3.4 


SM 


178 


Bao et al. 


0.9 




3.4 


MD 


179 


Zhou et al. 


0.8 


0.32 


3.4 


LDA 


180 



The early researches on the elasticity of carbon nan- 
otubes are focused on their Young's modulus Y and Pois- 
son ratio v. A SWNT is a single layer of carbon atoms. 
What is the thickness h of the atomic layer? It is a 
widely controversial question. Three typical values of 
the thickness listed in Table III are adopted or obtain ed in 
the previous literature Refs. Hfl EH El and[l6MIl The 
first one is about 0.7 A obtained from fitting the atomic 
scale model with the elastic shell theory of 3D isotropic 
materials. 167 ! 168 ! 169 ! 170 ! 171 ! 172 ! 173 The second one is about 
1.4 A derived from molecular dynamics or finite element 
method. 174 ' 175 The third one is about 3.4 A adopting the 
layer distance of bulk graphite. 176 ! 177 ! 178 ' 179 ' 180 Recently, 
Huang et al. have investigated the effective thickness of 
SWNTs and found it depends on the type of loadings 

The size- or chirality-dependent elastic properties of 
SWNTs have also been discussed by molecular mechan- 
ics mode l 182 ' 183 ' 184 and ab initio calculations! 173 ' 185 The 
common conclusion is that the Young's modulus and 
Poisson ratio depend weakly on the diameter and chi- 
rality of SWNTs if the diameter is larger than 1 nm. 
Only for very small SWNTs, the size and chirality effect 
is evident. The SWNTs synthesized in the laboratory 
have usually the diameters larger than 1 nm; thus the 
size and chirality effect can be neglected safely. 

The axial tension properties of MWNTs depend on the 
layer number of MWNTs for the small layer number and 
approach quickly to the properties similar to the bulk 
graphite^^d 87 

The buckling and stability of carbon nanotubes under 
pressure or bending is a hot topic in the recent researches, 
where the critical pressure, moment or the equivalent 
quantity, critical strain, are highly concerned. A long 



enough carbon nanotube under an axial loading might 
be regarded as a Euler rod and the axially critical strain 
is 3 - 

elf = air 2 l/AL 2 oc (p/L) 2 , (136) 

where L, p and A are the length, radius and cross- 
sectional area of the carbon nanotube, respectively. X 
is the moment of inertia of the nanotube. The value 
of a depends on the boundary conditions of the carbon 
nanotube. This relation has been investigated by atomic- 
scale finite element metho d 188 ' 189 ' 190 and molecular dy- 
namics method or ab initio calculations . 191 ' 192 ' 193 ' 194 ' 195 
The basic numerical result is that the tube exhibits 
rod-like buckling behavior as the right-handed side of 
Eq. (|136p if L ^> p. The Timoshenko beam theory, a 
more complicated theory than Euler rod theory, is also 
employed to discuss the buckling of MWNTsJ^ 6 - The dif- 
ference between the results of both theories vanishes for 
large value of L/p. 

For a short carbon nanotube under axial loading, the 
continuous shell model of 3D isotropic materials are 
widely used^d^^o^d^^^L^^l The axially 
critical strain of a short SWNT if& 

ef c M = [a/V3(l-^)](Vp) oc p~\ (137) 

where p and h are the radius and effective thickness 
of the SWNT, respectively, v is the Poisson ratio of 
the SWNT. The value of a depends on the bound- 
ary condition of the carbon nanotube. For a short 
MWNT, the above relation is applicable for the outmost 
layer of the tube because the inter-layer interaction of 
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MWNTs is very smallj 2 ^ . It has also been investigated by 
atomic-scale finite element method) 188 ' 189 ' 190 ' 198 molec- 
ular dynamics method , 40 ' 192 ' 193 ' 194 ' 199 and nanoindent 
experiment i 201 1 202 It is found that the tube displays in- 
deed the shell-like buckling behavior as the right-handed 
side of Eq. ()137|) for the tube aspect ratio L/p < 10. 

The stability of a long SWNT under radial hydrostatic 
pressure might also be described by the continuous shell 
model of 3D isotropic materials, and the critical pressure 

j : 138 

vt M « P~\ (138) 

where p is the radius of the SWNT. This relation has 
recently been confirmed by Hasegawa and Nishidate^S 3 . 
through ab initio calculations. The stability of a MWNT 
under radial hydrostatic pressure might also has the sim- 
ilar relation as Eq. (|138p if only we take p as the outmost 
radius of the MWNT, because the transverse elasticity of 
MWNT o 204 ' 205 is much weaker than the in-plane elastic- 
ity of the outmost single layer of tube. 

Bending can also result in the buckling of SWNTs. 
The kink phenomenon in a SWNT under pure bending 
has been investigated through molecular dynamics simu- 
lations and finite element methodj 40 ' 206 The critical cur- 
vature can be described as 

Kcr = e s z h c ell / P ocp- 2 , (139) 

where p is the radius of the SWNT. The kink phe- 
nomenon in a MWNT under pure bending satisfies the 
similar relation to Eq. Q139[) with small correction due to 
inter-layer van de Waals interaction o 207 ' 208 ' 209 if only we 
take p as the outmost radius of the MWNT. 

Here we would not further discuss the prob- 
lems on the buckling of MWNTs embedded in 
an elastic medium) 210 ' 211 ! 212 ! 213 ! 214 ! 215 the post- 
buckling behavior and the plastic properties 
of carbon nanotubes^^^ 2 ! 9 .* 2 ^ 2 ^ 2 . a s 
well as the mechanical properties of nanotube 
composites ) 223 ' 224 ' 225 ' 226 ' 227 ' 228 rather than recom- 
mend gentle readers to consult the corresponding 
literature. 



2. What are the fundamental quantities for SWNTs? 

As mentioned above, different thickness leads to dif- 
ferent Young's modulus (see Table |T]), which implies 
that the Young's modulus and thickness of SWNTs are 
not well-defined physical quantities! 2 ^ 9 - However, the in- 
plane Young's modulus Y s — Yh has the similar value 
22 eV/A 2 . Thus it is a more well-defined quantity than 
the Young's modulus and the thickness. Here we may 
ask: what are the fundamental quantities for SWNTs? 

A SWNT is also a single layer of graphite, whose de- 
formation energy can be also described as the revised 
Lenosky model (|121[) . The corresponding continuum 
limit is Eq. (|123p which contains four elastic constants 



k e , k, kd, and k. These four quantities avoid the con- 
troversial thickness of SWNTs. We suggest to use them 
as the fundamental quantities for SWNTs from which we 
can obtain some reduced quantities as follows. 

Let us consider a cylinder under an axial loading with 
line density / along the circumference. The correspond- 
ing axial and circumferential strains are denoted as en 
and £22 ■ With Eq. Q123|) , the free energy of this system 
is written as 

T « 2n P L[(k d /2)(e 11 + e 22 ) 2 - ten £22 - /e u ] (140) 

where L and p are the length and radius of the SWNT. 
The in-plane Young's modulus and Poisson ratio can 
be defined as Y s — //en and v s = —e^/sii- From 
dJ- '/en = and dT /e22 = 0, we derive 

Y s = k(2 - k/k d ) = 22.35 eV/A 2 , (141) 
v s = 1 - k/k d = 0.165, (142) 

where the value of Y s is close to the in-plane Young's 
modulus derived from Table |TJ It is in between 20- 
23 eV/A 2 obtained by Sanchez-Portal et al.l£l It is much 
larger than the value 15 eV/A 2 obtained by Arroyo et 
alML and Zhang et al.^- and 17 eV/A 2 by Caillerie et 
al.r^ but smaller than 34.6 eV/A 2 for armchair tube by 
Wang,— The value of v s is close to the value 0.16-0.19 
obtained by Yakobson et al.^ Kudin et ai.jiSI Pantano 
et a/. ; 170 ' 171 Wang et al.^^ Hernandez et ai.jiZS and Shen 
et aim. 

The other quantity, the bending rigidity D, is also 
widely discussed in literature. In terms of Eq. (|123p . the 
energy per area of a SWNT without the in-plane strains 
can be expressed as 

G g = k c /2p 2 = D/2p 2 . (143) 

Thus the bending rigidity 

D = k c = 1.62 eV, (144) 

which is quite close to the value 1.49-1.72 eV obtained 
by Kudin et al^- and Sanchez- Portal et al£3° through 
ab initio calculations. It is a little larger than the val- 
ues 0.85-1.22 eV obtained Yakobson et al.^ Pantano et 

aLr 170^71 and Wang ^34 

In terms of Eqs. (|141[) — (I144|) . we can infer the values 
of kd, k, k c from the previous literature, which are listed 
in Table [TTJ There is still lack of literature on k except 
our previous wor k 42 ' 44 ' 154 and the present review. More 
work on k would be highly appreciated in the future. 

We should emphasize that our formula (|123[) holds ap- 
proximate up to the order of (ro/p) 2 for SWNTs, where 
r is the C-C length and p the radius of the SWNT. The 
omitted terms is in the order of (ro/p) 4 . This is the main 
reason for the size effect on the el astic constants in the 
very small SWNTs found in Refs. fl7lll82Ml83lfl8ll85l . 
Additionally, we have not considered the effect of Stone- 
Wales defects on the loc al pr oper ties of carbon nan- 
otubes. In terms of Refs. |235| and |236| . we can deduce 
that the defects reduce the the elastic constants of car- 
bon nanotubes. 
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TABLE II: The values of Y s , v a , kd, k, k c and k. (MD = molecular dynamics; TB = tight-binding; SM = structure mechanics; 
FEM = finite element method; LDA = local density approach; CTIP=continuum theory of interatomic potential) 



Authors 


Y s (eV/A 2 ) 


v s 


k d (eV/A 2 ) 


k (eV/A 2 ) 


k c (eV) 


k (eV) 


Method 


Refs. 


Yakobson et al. 


22.69 


0.19 


23.54 


19.06 


0.85 


- 


MD 


40 


Tu & Ou-Yang 


22.03 


0.34 


24.88 


16.44 


1.17 


0.75 


LDA 


44 


Tu & Ou-Yang 


21.63 


0.18 


22.35 


18.33 


1.30 


0.88 


LDA 


154 


Kudin et al. 


21.69 


0.15 


22.19 


18.86 


1.49-1.53 


- 


ab initio 


167 


Zhou et al. 


23.59 


0.24 


25.03 


19.02 


1.14 


- 


TB 


168 


Pantano et al. 


22.5 


0.19 


23.34 


18.91 


1.09 


- 


SM & FEM 


170.171 


Chen and Cao 


34.38 


- 


- 


- 


- 


- 


SM 


172 


Wang et al. 


21.36 


0.16 


21.92 


18.41 


0.82 


- 


ab initio 


173 


Sears & Batra 


20.94 


0.21 


21.90 


17.30 


3.28 


- 


MD 


174 


Tserpes et al. 


22.05 


- 


- 


- 


- 


- 


FEM 


175 


Lu 


21.25 


0.28 


23.06 


16.60 


- 


- 


MD 


41 


Hernandez et al. 


25.50 


0.18 


26.35 


21.61 


- 


- 


TB 


176 


Shen & Li 


23.38 


0.16 


23.99 


20.16 


- 


- 


force-field 


177 


Li & Chou 


21.25 


- 


- 


- 


- 


- 


SM 


178 


Bao et al. 


19.13 


_ 


_ 


_ 


_ 


_ 


MD 


179 


Zhou et al. 


17.00 


0.32 


18.94 


12.88 






LDA 


180 


Sanchez-Portal 


19.41-22.40 


0.12-0.19 


19.92-23.00 


16.73-19.31 


1.49-1.72 




ab initio 


230 


Arroyo et al. 


15.19 


0.40 


18.08 


10.85 


0.69 




FEM 


231 


Zhang et al. 


14.75 












CTIP 


232 


Caillerie et al. 


17.31 


0.26 


18.57 


13.74 






CTIP 


233 


Wang 


34.63 or 17.31 








1.12 or 1.21 




CTIP 


234 


Present work 


22.35 


0.16 


22.97 


19.19 


1.62 


0.72 


LDA 





3. Revisit the stability of SWNTs 

Now we will revisit the stability of SWNTs with the 
four fundamental quantities k c , k, kd, and k or the cor- 
responding reduced quantities. 




FIG. 15: Bent SWNT. p and l/« are the radii of the SWNT 
and the centerline of the SWNT, respectively. 

First, let us consider a bent SWNT as shown in Fig. IT51 
where p and 1/k are the radii of the SWNT and the 
centerline of the SWNT, respectively. Assume that the 
centerline of the SWNT is not extended and the cross 
section of the SWNT is still flat after bending under the 
condition p <C L <C 1/k, where L is the total length of 
the centerline. In terms of Eq. (|123p . we can derive the 



deformation energy due to bending as 

A.F« f (k rod /2)n 2 ds, (145) 
Jo 

where ds is the arc length element of the centerline. 
The bending modulus of the rod is k ro d = 7tp[(2 — 
k/k d )kp 2 + k c ]. For the SWNT with diameter in the 
order of 1 nm, we can estimate k c <C (2 — k/k d )kp 2 . 
Considering Eq. (|14ip . we have 

k rod ^irY s p 3 . (146) 

If an axial compression force F is loaded on the both 
ends of the SWNT, following Euler rod theory^ we can 
easily derive the critical force, above which the SWNT is 
instable, as 

F r c od = 2airk rod /L 2 , (147) 

where a depends on the boundary conditions in two 
ends of the SWNT. Defining the critical strain as e r ° d = 
F£ od /2irpY s and considering the above two equations, we 
can derive 

e r z ° d = a(p/L) 2 . (148) 

This relation has the same asymptotic behavior as 
Eq. (|136p . which, as mentioned above, has been con- 
firmed by a lot of theoretical and numerical researches. 
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Secondly, let us consider a short SWNT with radius p 
and an axial compression force loaded on its two ends. 
The force per length along the circumference is denoted 
as /. Following Ru's workji^l considering Eq. (|123[) we 
have the critical axial force density as 



f! heU = a^Ws/p, (149) 

and the corresponding critical strain 

ef c M = f c /Y s = {a/p)^k~JY s , (150) 

where a depends on the boundary conditions in two ends 
of the SWNT. Y s is the in-plane Young's modulus as 
shown in Eq. (|14ip . The above relation (|150[) has the 
same asymptotic behavior as Eq. (|137p . which, as men- 
tioned in above, has been confirmed by several theoretical 
and numerical researches. 

Thirdly, let us consider a long enough SWNT with 
radius p and a radial compression pressure p loaded on its 
surface. In terms of the similar method on the stability 
of cell membranes, we can derive the critical pressure 

V t M = 3fc c /p 3 , (151) 

above which the SWNT will lose its stability. This rela- 
tion has the same asymptotic behavior as Eq. (|138p . The 
corresponding critical circumferential strain is 

ef = 2ppf r ell /Y s = 6k c /Y sP 2 . (152) 

Comparing Eq. (|151[) with ([55)1 . one can find that the 
critical pressures for carbon nanotubes and lipid tubules 
are in the same form. Yin et a l. have noticed this sim- 
ilarity in the recent work |237| . However, the profound 
mechanism is different: nanotubes can endure the shear 
strain while lipid tubules cannot. 

It seems that no literature discusses the possible in- 
stability of a SWNT under axial tension. Here we will 
give a qualitative analysis. Assume the tension density 
(i.e., force per length) to be /. The axial strain under 
the tension is f /Y s and the corresponding circumferen- 
tial strain is v s f /Y s . When it is beyond the critical value 
(fH)2"]) . the SWNT will be instable. Thus we obtain the 
critical tension density 

fr = 6k c /v sP 2 . (153) 

Only if f c sn is below the strength of the SWNT, the 
buckling phenomenon under tension can be observed. 

Till now, we have not found that k exists explicitly 
in the above equations (|148[) - (|152p for nanotubes. Be- 
cause the term related to k in the free energy (|123p can 
be transformed into the boundary term with the aid of 
Gauss-Bonnet formula, k should be implicitly contained 
by a in these equations, which need the further investi- 
gations in the future. 



V. CONCLUSION AND PROSPECT 

In summary, we present the elastic theory of low- 
dimensional (one- and two-dimensional) continua and its 
applications in bio- and nano-structures. The elastic the- 
ory of Kirchhoff rod, Helfrich rod, bending-softcn rod, 
fluid membrane, and solid shell is revisited. We construct 
the free energy density of the continua on the basis of 
the symmetry argument. The fundamental equations can 
be derived from the bottom-up and the top-down stand- 
points. Although they have different forms obtained from 
these two viewpoints, several examples reveal that they 
are, in fact, equivalent to each other. We investigate the 
kink stability of short DNA rings, the elasticity of lipid 
membranes, and the adhesions between a vesicle and a 
substrate or another vesicle. A cell membrane is simpli- 
fied as a composite shell of lipid bilayer and membrane 
skeleton. The membrane skeleton is shown to enhance 
highly the mechanical stability of cell membranes. We 
propose a revised Lenosky lattice model based on the lo- 
cal density approximation and derive its continuum form 
up to the second order terms of curvatures and strains, 
which is the same as the free energy of 2D solid shells. 
The intrinsic roughening of graphene and several typi- 
cal mechanical properties of carbon nanotubes are ad- 
dressed by using this continuum form. We can aban- 
don the controversial thickness and Young's modulus of 
graphene and SWNTs if we adopt this continuum form 
to describe the mechanical properties of graphene and 
SWNTs. 

Finally, we would like to list a few open problems which 
need to be addressed in the future work. 

(i) The vesicles with lipid domains have been investi- 
gated in Sec. IIII Bl There is a special lipid domain at 
liquid-ordered phase, so called the raft, which is enriched 
in cholesterol and sphingolipids. Cholesterol is a kind 
of chiral lipid molecules, which has not been included in 
the previous and present theory of lipid domains. A new 
theory with the chirality on the raft domain should be 
developed. 

(ii) The composite shell, as a model of cell mem- 
branes, has been investigated in Sec. IIII CI where the 
constraint between the area of the lipid bilayer and 
membrane skeleton is totally neglected. Additionally, 
only the small deformation of cell membranes are ad- 
dressed in this review. The large deformation behav- 
ior of cell membrane o 238 i 239 i 240 has recently discussed 
through numerical simulations. It is necessary to re- 
consider the composite shell model with the constraint 
J JdA = J Jj,dA and its behavior under large deforma- 
tion theoretically. 

(iii) We suggest adopting four parameters k c , k, k c i, 
and k to describe the mechanics of graphitic structures 
in Sec. IIVI However, there are sparse studies on k in 
previous literature. It is highly expected to theoretical 
and experimental work on this quantity. 

(iv) We only talk about the thermal fluctuation on the 
discussion of graphene. The fluctuations of DNA, lipid 
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membranes, and cell membranes are not in the range 
of our topics, on which we suggest that gentle readers 
consult Refsll3l3il63l63. 

(v) The elastic theory presented in this review is a 
static theory. Thus we are very regretted that we have 
to omit several important subjects such as the vesicles 
in shear flows , 241 i 242 i 243 i 244 i 245 an d dynamic response of 
carbon nanotubes or nanotube network o 246 ' 247 and so on. 
These topics will be quickly developed in the future. 
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Note added in proof 

After this review was in press, we noticed that re- 
cent researches 248,249 on the mechanical properties of 
nanosprings 248 and amorphous straight nanowires 249 
within the framework of Kirchhoff rod. We were also in- 
formed of the researches 250,251 by Arroyo and Belytschko 
on discussing the buckling pattern of multi- walled carbon 
nanotubes under pure bending. Additionally, we empha- 
size that the similar equations to (114) and (115) without 
p first obtained by Zhang et al. 252 



APPENDIX A: CURVE VARIATIONAL THEORY 

Here we sketch merely the derivation of Eqs. (|16p and 
(fT"7)) because the derivation of Eq. (fT5|) is trivial. 

For convenience, we denote ei = N, e2 = B, e3 = T. 
Then we have dr = 10363 and dei — LOij&j where 

u> 3 = ds, lo 12 = tw 3 , w 13 = -K0J3, u 2 3 = 0. (Al) 

Following the spirit of Ref. 14^ . any infinitesimal deforma- 
tion of a curve can be achieved by a displacement vector 
at each point on the curve as 

5r = v = fiiei + fi 2 e 2 + ri 3 e 3 , (A2) 

where 5 can be understood a variational operator. The 
frame is also changed because of the deformation of the 



curve, which is denoted as 

'•*'/ «->,,( V (i = 1,2,3), (A3) 

where f2y = — fiji, = 1,2,3) corresponds to the ro- 
tation of the frame due to the deformation of the curve. 
From Sdr — d8r, 6dej = dSej, and 6d<fr = d8(j> = 0, we 
can derive 

8UJ3 = + VL2UJ23 + dfl 3 , (A4) 

w 3 r2 3 i = dQi + fl 2 uj2i + Sl 3 w 3 i, (A5) 
w 3 fi 32 = + dfl 2 + Sl 3 w 3 2, (A6) 

8u)ij = dflij + flauij - ujuClij, (A7) 
5<j>' ds = -<t>'5u 3 . (A8) 

Considering the above equations (|A1[) . (|A4[) — (|A8|) and 
integral by parts as well as Stokes' theorem, we can derive 
Eqs. (HU) and (fTZj) from the free energy (fH|). 

APPENDIX B: SURFACE VARIATIONAL 
THEORY 

Following the spirit of Ref. 0, any infinitesimal defor- 
mation of a surface can be achieved by a displacement 
vector at each point on the surface as 

Sr = v = f^iei + fi2e2 + £! 3 e 3 , (Bl) 

where 5 can be understood a variational operator. The 
frame is also changed because of the deformation of the 
surface, which is still denoted as 

Se, <>,,<>,. (1 = 1,2,3), (B2) 

where fiy = — flji, (i, j = 1,2,3) corresponds to the ro- 
tation of the frame due to the deformation of the surface. 
From Sdr = dSr, Sdej — dSej, and Eqs. ©-(O, we can 
derive 



5uj\ 


= dv ■ ei - W2O21, 


(B3) 


Su>2 


= d~v ■ e2 — W1O12, 


(B4) 


r>l 3 


= £l 3j i + afli + bil 2 , 


(B5) 


^23 


= Q 3t 2 +bQ 1 +cQ 2 , 


(B6) 


StUij 


= dClij + Oj/OJy — WyOy. 


(B7) 



These equations are the essential equations of the surface 
variational theory based on the moving frame method. 
With them as well as Eqs. {5} and ([6]), we can easily 
derive 

SdA = (div v - 2HQ 3 )dA, (B8) 
5{2H) = [V 2 + {AH 2 - 2K)]Q 3 + V(2H) ■ v, (B9) 
SK = V ■ VQ 3 + 2KHn 3 + VK ■ v. (B10) 

Using the above three equations and the Stokes' theorem, 
we can easily derive Eq. (|57l) from the free energy (|56|). 
or Eqs. (|87 p - (150"l) from the free energy ([55]) . and so on. 
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APPENDIX C: STOKES' THEOREM AND THE 
OTHER IMPORTANT GEOMETRIC RELATIONS 

The Stokes' theorem is a crucial theorem in differential 
geometry. Let us denote the boundary of domain D as 
81). The Stokes' theorem states as: If oj is a differential 
form on dD. then 



du. 



(CI) 



In particular, J s du> = for a closed domain D. 

It contains a lot of geometric relations, which are listed 
as follows. 

(i) For any smooth functions / and h on 2D domain 
2), we have 



f {fd*dh- 


- hd*df) = 


I (f*dh- 
Jav 


h*df), 


(C2) 


I (fd*dh- 


hd*df) = 


i (f*dh- 


h*df), 


(C3) 


/ (fd*dh- 


- hdUf) = 


I (fMh- 


h*df). 


(C4) 



where d and * are generalized differential operator and 
Hodge star which satisfy df = fiUis + /2W23 and *df = 
/1W23 - /2W13 if df = fxui + / 2 w 2 -— 

(ii) If u is a vector defined on a closed surface, then 

du ■ A * du = — J u • d* du, (C5) 

where the dot represents the inner product of vectors. 



(hi) For the tensors 6 and 9Jt defined in Sec. [HI we 
have 



6 • b ds = / div 6 dA, 
Wl-hds = / divVJldA, 



(C6) 
(C7) 



where J) is a 2D domain with boundary dD and b is the 
normal vector of dD in the tangent plane. 

The above three items are also called the Stokes' the- 
orem in this review, which are widely used in the varia- 
tional process. 

The other geometric identities linking the vector form 
and differential form on a smooth surface used in this re- 
view are summarized as follows without additional proof. 



(curl u) dA = d{u ■ dr), 


(C8) 


(div u) dA — d(*u ■ dr), 


(C9) 


(V ■ u) dA = d(*u ■ dr), 


(CIO) 


(V ■ u) dA = d(*u ■ dr), 


(Cll) 


Vf-dr = df, 


(C12) 


Vf-dr = df, 


(C13) 


(V 2 /) dA = d * df, 


(C14) 


(V ■ V/) dA = d * df, 


(C15) 


(V • V/) dA = dUf, 


(C16) 


(V/ • u) dA = df A *u • dr, 


(C17) 


(V 2 u) dA = d * du, 


(C18) 


(Vu) • dr = du. 


(C19) 
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